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ABSTRACT 



PHYSICS OF SINGLE-TOP QUARK PRODUCTION AT HADRON COLLIDERS 

By 

Douglas Olaf Carlson 

We discuss the physics of single-top quark production and decay at hadron col- 
Uders, such as the Tevatron, the Di-TeV and the LHC. Our study includes how to 
measure the mass and the width of the top quark produced from a single-t or a single-f 
process. We also show how to probe new physics by studying the couplings of t-b-W 
and show what can be improved from measuring the production rate of single-top 
events. We also discuss how to probe CP properties of the top quark by measur- 
ing the single-top production rate. Finally, we present a Monte Carlo study on the 
detection of single-top events in hadron collisions. 
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iii 



ACKNOWLEDGEMENTS 



I would like to express my deepest gratitude to the people who have contributed 
to me and to this work: 

To my thesis advisor, C.-P. Yuan, for his calm guidance over the past three years; 

To my mentor. Chip Brock, for sparking my interest in this field and introducing 
me to C.-P.; 

To the members of my Thesis Committee: Chip Brock, Wu-Ki Tung, S. D. Ma- 
hanti and Horace Smith, for their carefuU reading of the manuscript; 

To my colleagues in the High Energy Theory group: Glenn Ladinsky and Pankaj 
Agrawal for many illuminating discussions; and my graduate student colleagues: Ehab 
Malkawi, Mike Wiest, Csaba Balasz, Liang-Hung Lai and Xiaoning Wong for our 
lunch gatherings; 

To my wife. Dawn, for her love and support; 

To my parents, Jerry and Justine, for everything. 



iv 



Contents 



LIST OF TABLES 
LIST OF FIGURES 

1 Introduction to the Standard Model 

1.1 The Matter Spectrum of Spin-^ Fermions 

1.2 Force Mediators as Spin-1 Gauge Bosons 

1.3 Gauge Invariant Interactions of Fermion and Gauge Boson Fields 

1.4 Gauge Field Self Interactions" 

1.5 'I'he Higgs Mechanism 

1.6 Review 

2 Introduction to the Top Quark 

3 The Single- Top Production Mechanism 

4 Measuring the Top Quark Mass 



5 Measuring the Top Quark Width 



6 Top Quark Couplings to the W Gauge Boson 



3.1 i'rom the Decay of 'Ibp Quarks 



3.2 From the Production of 'Ibp Quarks 



7 Probing CP Properties in Top Quarks 

8 A Monte Carlo Study 





'ievatron with a 




= 2TeV 




B.2 


'ievatron with y 




= 4TeV 



p — LHC with V^= l4TeV 



9 Discussions and Conclusions 



A Helicity Amplitude Method 



Helicity Amplitudes for ub dt 



A. 2 Helicity Amplitudes for ub —>■ dt in the CMS 



A. 3 Helicity Amplitudes for ug — > dt(- 


-^bW+i- 




b 










A.4 Helicity Amplitudes for M 6 ^ (it(- 


-^bw-^i- 







A. 5 Helicity Amplitudes for ud W* ^ btj^ bW+j^ i+ Uj)) 
A. 6 Helicity Amplitudes for ud ^ fefePF "*"(—> v^) 



B Event Rate of the (2 — > 3) Process ug — > dtb 



C Helicity Amplitudes of t ^ W^b and t b 
D The Total Rate for l4/-gluon l^Yision 



vi 



E The Eikonal Approximation for a{AB{ud) — » bbW^ + jet) 
1^' The Computer Program ONETOF" 



LIST OF REFERENCES 



vii 



List of Tables 



1.1 Lepton and Quark Masses 

1.2 Boson Masses 

1.3 Quantum numbers of the fermion spectrum 

1.4 Quantum numbers of the Higgs doublet 

3.1 Rates of the above processes for rrit = 180(140) GeV. (Branching ratios 

are not included here.) For VS* = 2TeV and 4TeV we include rates 

for a pp machine. At y/S = 14 TeV the rates are for a pp machine! 
For the single-top rates we only include single-t production. 



3.1 Results on the accuracy of measuring f^'^ for various luminosities! 
(Only statistical errors are included at the 95% confidence level.) . . 



viii 



List of Figures 



2.1 Diagrams contributing to the QCD production of qq, gg — » tt 



3.1 Diagrams for various single-top quark processes 



3.2 Rate in 


pb 


for qq, gg tt, q'g{W'^g) — > q 


tb, q'q W* tb and 


gb ^ W t at various energies of pp colliders. 




3.3 Rate in 


pb 


for qq, gg tt, q'g{W^g) qtb, q'q ^ W* ^ tb and 


gb ^ W t at various energies of pp colliders 



11 

14 
15 
16 

18 

19 

20 

21 

25 
25 

27 

33 
33 



3.4 I'eynman diagrams illustrating the subtraction procedure for calcu- 



lating tiie total rate -gluon fusion: q'b ^ qt (B q'g{W~*^ 



qtb e {a 




bb ® q'b ^ qt) 




3.5 Rate in 


pb 


for single-t production: q'b qt (2 — » 2), q' g — > qtb 


V^' 


and the splitting piece g bb ® q'b ^ qt in which bb are 


coUmear. 


'I'he rates are for pp colliders 



3.6 Rate in 


pb 


for single-t production: q'b qt 


(2 - 2), q'q 


qtb 


(2^3) 


and the splitting piece g ^ bb ® q'b - 


-* qt in which bb 


are 


collinear. 


The rates are for pp colliders 







3.7 Rate of PF-gluon fusion process versus scale Q for = 180 GeV and 
V^ = 2TeV 



1.1 I'he lepton-|-jet decay mode of tt production. 



O I'he di-lepton decay mode of tt production. 



4.3 Distributions of mu (solid) and mj,f (dash) in tt events for a 180 GeV 
top quark. 



|6.1 For a left-handed t-b-W vertex. 



|6.2 For a right-handed t-b-W vertex. 



3.3 nihi distribution for SM top quark (solid) and for pure right-hande^ t-b-W 
coupling of tbW (dash). 34 

3.4 cos Of distribution for SM top quark (solid) and for pure right-hand) t-b-W 
coupling of t6 14/ (dash). 35 



ix 



3.5 Constraint on 




and K/j*"^ given 


20% and 50% error in measure- 


ment of Standard Model rate for W-g 


uon fusion. Curves are identical 


for mt = 140 GeV and nit = 180 GeV 



g.l Diagrams for ud, cs — » bbW'^{—>- E^uJT] 45 



j.2 'I'he rapidity distribution of the spectator quark q, after cuts in Equa"^ 
+i — I a c — — „i , „-*-/" , , /)+,,\\ „„J „f i-u„ „ — „ 



tion (8.3), foj: the signal q'b — > > bW~^{^ of the spec- 
tator quark b for the major background q'q —>■ bbW~^{^ i'^u) (dots), 
for nif = 180 GeV (sohd) and 140 GeV (dash), at the Tevatron. (The 
vertical scale is arbitrary, but the relative size among these curves are 
absolute.) 46 



B.3 Pt distribution of the b quark, after requiring > 15 GeV along with] 



all the other cuts in (8.3| ), for the signal q'b — » qt{^ bW^{—^ i^J^)), 



and the major background q'q — > bbW^{^ at the Tevatron. . . 48 

g.4 The cosOjc, distribution prior to cut Equation (8.4| ) for the signal q'b —>■ 



qt{-^ bW^{-^ d-^y)) and the major background q'q bbW^{-^ d~^^), 

at the Tevatron 49 



g.5 The mt distribution after the cuts Equation (8.3| ) and Equation (^]4| 



for mt = 180 GeV (solid) and 140 GeV (dash) at the Tevatron including 
both the signal and background events with — >■ or /i^ 50 



g.6 Same as Figure 8.5|, but with detector resolution effects as described 



in Equation 3.9 53 



B.7 mu distributions without (top) or with (bottom) smearing, after all 



the above analysis for the Tevatron. 54 



3.8 The rapidity distribution, after Equation (8.3|), of the spectator jet 



[i.e., 6-jet if 6-jet identified) in the W* event for mt=180 GeV (solid) 

and 140 GeV (dash) against Figure ^]2] (dots), at the Tevatron. ... 55 



3.9 'I'he production rate for a left-handed (long dash) or a right-hande3 



(short dash) top quark from tlie W* process. TUT upper solid line Is] 
tET total rate for the W-gluon fusion process, the lower solid line foi] 



process.^ 60 



3.10 'i'he distribution of cosOig in W* 


event for mt=lSO GeV (solid) and 


140 GeV (dash) against Figure 8.4 


(dots), at the Tevatron 



3.11 ARhh distributions in W + bb + jet (solid) and W + bb (dash) events 



after applying the cuts listed in (8.11| ) 62 



3.12 Mf,f, distributions in W + bb + jet (solid) and W + bb (dash) events 



after applying the cuts listed in (8.11| ). 63 



3.13 The rapidity distribution of the spectator quark q, after cuts in Equa- 
tion (8.12|) , for the signal q'b — » qtj-^ bW'^i—^ d'^i^)), and of the spec- 



tator quark b for the major background q'q —> bbW^{^ i^u) (dots), 
for mt = 180 GeV (solid) and 140 GeV (dash), at the Di-TeV 64 



3.14 Pt distribution of the b quark, after requiring P^ > 15 GeV along with] 



all the other cuts in (8.12| ), for the signal q'b — > qt{^ bW^{-^ ^^^)); 



and the major background q'q — * bbW^{—>- l^v), at the Di-TeV. ... 65 



X 



g.l5 cosOfg distribution for the signal q'b — > » bW'^{—^ ^^^)) aiid back" 



ground q'q —>■ bbW 



z/) at the Di-TeV. 



g.l6 The nit distribution after the cuts Equation (8.12|) and Equation ( ^.41) 
for = 180 GeV (sohd) and 140 GeV (dash) at the Di-TeV, including 
both the signal and background events with 



or fi"^. 



g.l7 Same as Figure 8.16| but with detector resolution effects as described 
in Equation ( |8.9| )" 

g.l8 Pt distribution of the b quark, after cuts in (8.13| ) , for the signal q'b 
qt{-^ bW^{-^ i^T^)), and the major background q'q — > bbW^ 
at the LHC 



V] 



j.l9 'I'he rapidity distribution of the spectator quark g, after cuts in Equa- 
tion (8.13| ), for the signal q'b — » > bW^{-^ of the spec- 

(dots) 



66 

67 
68 

69 

70 
71 

72 

73 

78 
83 
84 
85 
86 

87 

95 



tator quark b for the major background q'q —> bbW'^{—^ t 
for mt = 180 GeV (solid) and 140 GeV (dash), at the LHC. 



B.20 cosOia distribution for the signal q'b 
ground q'q 



EL 

a+u) at the Lmr. 



bW+{^i+iy)) and back- 



j.21 I'he rrit distribution after the cuts Equation (8.13| ) and Equation 

for mt = 180 GeV (solid) and 140 GeV (dash) at the LHC includmg 
both the signal and background events with 



or fi"^. 



g.22 Same as Figure 8.21| , but with detector resolution effects as described 
in Equation (8^) 



A.l Diagram for the (2 — » 2) process ub ^ dt 



A. 2 Diagrams for uq dt{^ bW~^{-^ d'^ ue))b. 



A. 3 Diagram for ^ 



A.4 Diagram for ud bt(^ bW^{^ i'^ i^e)) 



A.5 Diagrams for ud bb W^{^ Vj). 



B.l Diagrams for u g —>■ dtb . 



D.l n-body scattering 



xi 



Chapter 1 

Introduction to the Standard 
Model 



The Standard Model (SM) of elementary particle physics ||T], ||, ^ is a Yang-Mills 
gauge field theory with symmetry 

SU(3)c X SU(2)^ X U(l)y. (1.1) 

It has been very successful in explaining and predicting experimental data. The 
SU(3)(-^ sector governs the strong force of the SM and is known as quantum chromo- 
dynamics (QCD). The SU(2)^ x U(l)y sector unifies the electromagnetic and weak 
forces, collectively known as the electroweak force. The unification of the electromag- 
netic and weak forces is accomplished in the SM via the mixing of the neutral SU(2)^ 
gauge boson and the hypercharge gauge boson of U(l)y. Masses are introduced in 
a gauge invariant way through spontaneous symmetry breaking which gives rise to 
the as yet undiscovered Higgs boson (H). This process of spontaneous symmetry 
breaking is known as the Higgs mechanism. 

The standard model does not incorporate the gravitational force. So far, no 
definitive quantum field theory of gravity exists. Gravity affects all massive particles, 
however, gravitational interactions with elementary particles are too weak and can 
be ignored. 

In the standard model there are three generations of leptons and quarks as listed 
in Table LT| . Associated with each force is one or several gauge bosons as listed in 
Table |L2. For reference we list the masses of each particle as found in Reference 0. 
So far, the only particle for which a discovery is lacking is the scalar Higgs boson. 

In what follows, we briefly describe the particle spectrum, particle interactions 
and the 18 independent parameters which constitute the standard model. 



1.1 The Matter Spectrum of Spin-| Fermions 



The matter spectrum consists of twelve fermions which are organized in Table 1.1 



according to the symmetry structure of Equation ( |l.lj ). Each successive generation 



1 



2 



Table 1.1: Lepton and Quark Masses 



Particle 


Symbol 


Mass (GeV) 




Electron neutrino 









Electron 


e 


0.00051 


First 


Up quark 


u 


0.002 to 0.008 


Generation 


Down quark 


d 


0.005 to 0.015 




Muon neutrino 









Muon 




0.106 


Second 


Charm quark 


C 


1.3 to 1.7 


Generation 


Strange quark 


s 


0.1 to 0.3 




Tau neutrino 









Tau 


T 


1.78 


Third 


Top quark 


t 


174 


Generation 


Bottom quark 


b 


4.7 to 5.3 





Table 1.2: Boson Masses 



Particle Symbol Mass (GeV) 



Photon 
W Boson 
Z Boson 
Gluon 



7 
G 



Electromagnetic Force 

80.22 Charged Weak Force 

91.187 Neutral Weak Force 

Strong Force 



Higgs 



H 



60 < niH < 800 Spontaneous Symmetry Breaking 
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Table 1.3: Quantum numbers of the fermion spectrum 



Cliirality 


Q 




Y 




C 









1/2 


-1 









gl 


-1 


-1/2 


-1 









Ul 


2/3 


1/2 


1/3 


r, 


^, 




dL 


-1/3 


-1/2 


1/3 


r, 


9, 


b 




-1 





-2 









Ur 


2/3 





4/3 


r, 


^, 


b 


dR 


-1/3 





-2/3 


r, 


^, 


b 



is a more massive copy of the previous generation, so only the quantum numbers 
of the first generation are shown in Table |1.3| . Under the Sl]{2)^ sector of the SM, 
left-handed fermions transform as weak isospin (Tw) doublets, 




whereas right-handed fermions transform as singlets, 

(e)i?, iu)R, {d)R. (1.3) 

Since neutrinos are massless Dirac fermions in the SM, there are no right-handed 
neutrinos. Once the third component of Weak isospin T-^ is assigned, the values 
of hypercharge Y can be determined to cancel the chiral anomalies 0. With these 
quantum numbers in place, the charge quantization is determined by 

Q = T^ + j. (1.4) 



Only quarks transform under the SU(3)(^ sector of the SM. Each quark flavor, 
{i.e., q = u,d, s,c,b,t), carries red{r) , green{g) or blue{b) color charge and therefore 
transforms triplet: 




1.5) 
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1.2 Force Mediators as Spin-1 Gauge Bosons 

Associated with the group structure of the SM are twelve generators and each gen- 
erator is associated with a gauge boson. Therefore, U(l)y has one generator, the 
neutral hypercharge gauge boson 5^. SU(2)^ has three (2x2 — 1 = 3) generators, 
two of which are charged S\J{2)^ gauge bosons VF^ and one neutral SU(2)^ gauge 

boson W^. Finally, SU(3)^ has eight (3x3 — 1 = 8) generators and therefore, eight 
gluons, G;j,a = 1,2, ...,8. 

As stated earlier, electroweak unification is accomplished by mixing the 5^ and 
gauge bosons. Formally this is accomplished via the following rotation, 

where 6w, called the weak mixing angle, is chosen such that only couples with 
charged particles. is then identified as the photon field in quantum electrodynam- 
ics (QED) and additionally, a neutral weak force, Z° , is obtained. 

At this point in the theory, all fermions and gauge bosons are massless to preserve 
gauge invariance. In the next two sections we describe how the force mediators inter- 
act with fermions and amongst themselves. The Higgs mechanism is also introduced 
to incorporate mass in a gauge invariant way. 

1.3 Gauge Invariant Interactions of Fermion and 
Gauge Boson Fields 

We begin with the lagrangian for a massless free fermion field 

£^^ = -^^^. (1.7) 

Equation ( |1.7| ) is the kinetic term for fermion fields.Q To make the lagrangian gauge 
invariant, we introduce the gauge covariant derivative 

d,^D^ = d,- tg.^B, - tg2^W/, - tg-s^C;,, (1.8) 

where j = 1,2,3 and a = 1,2,..., 8. The r-^'s are the Pauli matrices defined in 

Appendix A and = ^ is the Weak isospin. The (3 x 3) matrices (generators) 

A", a = 1, 2, . . . , 8, are the SU(3)(^ matrices. As a result, we obtain gauge interaction 
terms in the lagrangian. 

The term acts on all fields with different y's, including leptons and quarks in 
Table (same for the other two generations) and the Higgs doublet field $ discussed 

^ ^ — 9^7^ where 9^ is the Lorentz invariant space-time derivative and 7'^ are the Dirac matrices. 
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in Section 1.5 below. The term acts only on the SU(2)^ doublets with non-zero 

T-yj and the field $. In the process of obtaining electroweak unification the coupling 
constants, gi and g2 become related through the weak mixing angle 6w via 

^ = tangly, (1.9) 
91 

where gi sin 9]y = e. The value of e is related to the fine structure constant a by 
a = e^/An. 

For electroweak interactions, the colors of the quarks have to be the same since 
SU(2)j^ X U(l)y does not act on color space. Also, the and fields do not 

induce quark or lepton fiavor changing at the Born level, {i.e. no u ^ c ^ t, 
d ^ s ^ b, Ue ^ Un Ur or e ^ fi ^ T transitions occur). However, charged 
current transitions via W^j^ do occur and for quarks, are not restricted by generation 

{i.e. u,c,t ^ d,s,b transitions are allowed). Lepton fiavor changes are restricted 
to generation due to massless neutrinos. Finally, the SM does not allow for direct 
lepton-quark transitions. 

The term in Equation ( |1.8D only acts on colored fermions, i.e. quarks. The 
coupling strength g^ = gs is universal for all colored quarks. Analogous to the 
fine structure constant a in electromagnetic theory is the strong coupling constant 
as = gs/4:iT in QCD. Here again, no fiavor changes occur. 



1.4 Gauge Field Self Interactions 

To complete the lagrangian for massless fermions and gauge bosons and their inter- 
actions, we must introduce the kinetic term for gauge bosons. 



Cgk = -\b,.b^' - - \q%Q'''\ (1.10) 



where 



B^, = d^B, - d,B^, (1.11) 
= d,Wi - d^Wi + g^e'^'WiW^, (1.12) 



Ql. = d,Gt - d,G; + gsr'^GlGl. (1.13) 

The Lie group structure constants e*-'^ and f"'^'^ are defined through the following 
anti-commutation relations, 



2"'y 



^e^^'y, ^,J,fc = 1,2,3 (1.14) 



6 



and 

= ifabc—, a,6,c = 1,2, ... ,8 (1.15) 

for SU(2) and SU(3), respectively. 

We note that the pure Yang-Mills terms, — |W*y>V*^^ and —jG^^Q"'^'^ contain 

factors that are trilinear and quadrilinear in WJ^ and G^. These Yang-Mills terms 
expand out partially as 

92 t'^^ {d^WDW^^'W^'' - ^ e'^^ e^'"^ WlW^W^^W^"" (1.16) 

and 

9s r''" {d^GDG^^G^" - ^ P""" f"^' GlGlC^^G^" (1.17) 

respectively, and correspond to self-couplings of non-abelian gauge fields. This is 
fundamentally different than in the abelian case where, in QED, photons do not 
directly couple with photons. It also accounts for the short distance interaction of 
the strong force, despite the fact that the gluon is massless. As can be seen in 
Chapter 2, the triple gluon interaction contributes to ti production via gg ti at 
hadron colliders. 



A" A^ 



1.5 The Higgs Mechanism 

The goal of the Higgs mechanism is to introduce mass to the particles in the SM in a 
gauge invariant way. We begin by defining a complex doublet scalar field $ composed 
of four real scalar fields if, ch^, 6^ and 6"^ where 



The quantum numbers of the Higgs field shown in Table 11.41 



The fields 0° and (p"^ = (0^ =F i(p'^)/V^ are the unphysical would-be Goldstone 
bosons associated with spontaneous symmetry breaking. They give rise to the masses 
of the gauge bosons and Z^. One physical field thus remains, which is the Higgs 
field H. The constant v ~ 246 GeV is the scale characterizing the symmetry breaking 
scale and is called the vacuum expectation value of $, where 



(<|.)o^(0|<l>|0)= ( f ). (1.19) 
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Table 1.4: Quantum numbers of the Higgs doublet 

Q Y C 

^^±H^ 1/2 -1 
-1 -1/2 -1 



The lagrangian for the Higgs sector is 



= (D^$)t(L)M$) _ ^($t$)2 _ ^2($t$), (1.20) 



with 

[ Y tO ^ ^ ^ 

$ 1.21 



Y 

d^-igi-B^~ig2—Wl 



as in Equation (|1.8| ) without the gluon interaction. 

If ;U < and A > then the minimum of the potential energy occurs at 



The Goldstone bosons 0^ and (jP are "eaten" by the vector bosons and Z^, 
respectively, where 

Mw = ):92V and Mz = (1.23) 
2 cos Uw 

Therefore, and have three polarization states; two transverse and one longi- 
tudinal. The massless photon and gluon have only the two transverse polarization 
states. The Higgs mass is given by 



yX. (1.24) 



To introduce fermion mass in a gauge invariant way one introduces the Higgs 
mechanism through Yukawa coupling interactions. For the first generation 

\/2m„ ^ V2md , . . . , 

^Yukawa = [Ul ^l) $ MiJ H (^L ^l) {-ir2^ ) dn 

V V 

■\f2me 

H [pi^ Cl) (— ir2$*) Cr + hermitian conjugate (1-25) 
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where m„, md and rrie are up quark, down quark and electron masses, respectively 
(neutrinos are massless) and 

-ir2$*=(^° "q^ ) = I IH^M! I . (1.26) 



As mentioned in Section 1.3, neutral currents coupled to 7, and G do not 
change flavor, although G changes color, but charged currents coupled to do 
change flavor. For leptons, the flavor change does not exceed generational bounds, 
due to the massless neutrino. However, there is a chance that an up quark, for 
instance, can change to a down quark, a strange quark or even a bottom quark. 
This is called quark mixing, which is due to the weak eigenstates (indicated by the 
subscript "Weak") of quarks being different than the mass eigenstates (indicated by 
the subscript "Mass"). 

By convention, the three charge 2/3 quarks u, c and t are unmixed: 







(u 




^ Weak 


C 









Mass 



1.27) 



All the mixing is therefore expressed in terms of a (3 x 3) unitary matrix V operating 
on the charge (-1/3) quarks d, s and b: 




I-AV2 A AX^{p-ir])\ fd\ 

-X 1 -AV2 s . (1.28) 

AXHl-p-^v) -AX' 1 J V^Mass 



The matrix V is known as the Cabibbo-Kobayashi-Maskawa matrix (CKM), which 
consists of three mixing angles and one phase. The second parameterization is due 
to Wolfenstein 0, where A ~ 0.22, A ~ 1, r/ ~ 0.5 and -0.4 < p < 0.2. CP 
violation, which is the violation of combined charge conjugation C and the parity 
transformation P, is characterized by the CP violating phase in p — irj. 



1.6 Review 

In this introduction to the SM, we have outlined how the SM is constructed based 
on gauge invariance and the Higgs mechanism. Although the SM has been successful 
in describing experimental data, there are 18 independent parameters which must be 
determined experimentally. These parameters are: 



9 

• nine fermion masses: mejmi^jmrjrnujmdjrncjrnsjmtjmb (neutrinos are mass- 
less) , 

• four CKM parameters: A, A, rj, p, 

• four electroweak parameters: e, ^^y, M^y, mjy, 

• the strong coupling constant: as- 



Chapter 2 

Introduction to the Top Quark 



If the SU(2) structure of the Standard Model (SM) holds, the top quark {t) has to 
exist as the weak isospin partner of the bottom quark (6) 0. If the coupling of 
t-t-Z is as predicted in the SM, then from LEP and SLAG experiments the mass 
of the top quark (mj) has to be larger than half of the Z-boson mass (~ 45GeV) 
independent of how the top quark decays. If the coupling of t-b-W is as predicted in 
the SM, then from the measurement of the total width of the ly-boson, by measuring 
the ratio of the event rates of pp — *• W{—> ip) to pp — > Z{—>- [|, the mass of 

the top quark (mt) has to be larger than 62 GeV independent of how the top quark 
decays From examining the radiative corrections to low energy observables, such 
as the p parameter^ which is proportional to at the one loop level [1^, rrit has 



to be less than about 200 GeV. Based upon analysis of a broad range of Electroweak 
data, the mass of the SM top quark is expected to be in the vicinity of 150 to 200 GeV 
11, 12, 13, 14 1 . Independently, from the direct search at the Tevatron, the top quark 



has been discovered and found to have mass of = 176 ± 8 (stat.) ± 10(sys.) GeV 
from CDF data [H, and mj = 1991^? (stat.) ± 22(sys.) GeV from D0 data 



For a heavy top quark, rrit is of the order of the electroweak symmetry breaking 

— 1/2 

scale V = {\/2Gf) = 246 GeV. In fact, recall the Yukawa coupling interaction, 
this time for third generation quarks, 

Cvukawa = {iLbL)<^tn + (hh) (-^T^**) 6«. (2.1) 

V V 

We see that 



V 



1, (2.2) 



proton is denoted by p and an anti-proton by p 
^ P ~ JPtS^W^' where Mw (or Mz) is the mass of W"^ (or Z) boson. 6w is the weak mixing 

angle, p has been measured to the accuracy of about 0.1%. 
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11 





Figure 2.1: Diagrams contributing to the QCD production of gg, gg tt 



for rrit = 175 GeV. Because the generation of fermion mass can be closely related 
to the electroweak symmetry breaking [|r^, |TB|, effects from new physics should be 
more apparent in the top quark sector than any other light sector of the electroweak 
theory. Thus, the top quark system may be used to probe the symmetry breaking 
sector. A few examples were discussed in Ref. [|l9l to illustrate that different models of 



electroweak symmetry breaking mechanism will induce different interactions among 
the top quark and the W- and Z-bosons. Therefore, hopefully through studying the 
top quark system one may eventually learn about the symmetry breaking sector of 
the electroweak theory. 

The most important consequence of a heavy top quark is that to a good approxi- 
mation it decays as a free quark because its lifetime is short and it does not have time 
to bind with light quarks before it decays [^. Furthermore, because the heavy top 



quark has the weak two-body decay t — »• bW~^, it will analyze its own polarization. 
Thus we can use the polarization properties of the top quark as additional observables 
to test the SM and to probe new physics. In the SM, the heavy top quark produced 
from the usual QCD process, at the Born level, is unpolarized. However, top quarks 
will have longitudinal polarization if weak effects are present in their production |pi|| . 
For instance, the top quark produced from the H^-gluon fusion process is left-hand 
polarized. With a large number of top quark events, it will be possible to test the 
polarization effects of top quarks. 

How to detect a SM top quark pair produced via the QCD processes qq, gg tt, 
as shown in Figure |2.1| , has been extensively studied in the literature |2^. In this 



paper we will concentrate on how to detect and^ study the top quark produced from 
the single-top quark processes q'g(W~^g) — > qtb, q'b — > qt, gb W~t, and q'q 
W* tb. For the single-top production we will only consider the decay mode of 
t — > bW^ bi'^u, with = e^orz/"*". (The branching ratio for this decay mode is 
Br = i.) 

The rest of this paper is organized as follows. In Chapter 3 we discuss the pro- 
duction rates of top quarks at hadron colliders. Following that, we will discuss in 
Chapters 4 and 5, respectively, how to measure the mass and the width of the top 
quark. In Chapter 6 we discuss what we have learned about the couplings of the top 
quark to the weak gauge bosons and show what can be improved from measuring 
the production rate of single-top quark events. We will also discuss in Chapter 7 the 
potential of the Tevatron as a pp collider to probe CP properties of the top quark 
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by simply measuring the single-top quark production rate. Finally, in Chapter 8 we 
present a Monte Carlo study on the detection of a single-top quark event in hadron 
collisions. Various unique features of the kinematics of the single-top quark signal 
will be discussed. Chapter 9 contains our conclusions. Based upon the results of the 
FNAL CDF and D0 , the mass of the SM top quark rrit is not likely to be lighter 
than, say, 140 GeV. Throughout this paper we will use rrit = 180 GeV (or 140 GeV) 
as an example of a heavy (or a light) top quark for our studies. 



Chapter 3 



The Single- Top Production 
Mechanism 



In this chapter we consider the production rate of a single-top quark at the Tevatron, 
the Di-TeV (the upgraded Tevatron) and the LHC (Large Hadron Colhder) colhders. 
In referring to single-top production, unless stated otherwise, we will concentrate 
only on the positive charge mode {i.e., only including single-t, but not single-f). 
The colliders we consider are the Tevatron (a pp collider) with the Main Injector at 

= 2 TeV, the Di-TeV (a pp collider) at 4 TeV and the LHC (a pp collider) at = 

14 TeV with an integrated luminosity of lfb~^, 10fb~^, and 100 fb~^, respectively.^ 

A single-top quark signal can be produced from either the W^-gluon fusion process 
q'g(W'^g) qtb (or q'b qt) p3|,j2l, the Drell-Yan type process q'q ^ W* ^ th 
(also known as 'W*" production) ||25|| , or Wt production via gh W~t ||26|| . The 
corresponding Feynman diagrams for these processes are shown in Figure |3.1| . 

In Figures and |3.3| we show the total cross sections of these processes for 
the Tevatron, the Di-TeV and the LHC energies referred to above. For reference 
we include plots of the cross sections of top quarks as a function of mt in both 
the pp collisions, shown in Figure |3.2|, and pp collisions, shown in Figure |3.3|. The 



parton distribution function (PDF) used in our calculation is the leading order set 
CTEQ2L [g^. We note that taking the Aqcd value given in CTEQ2L PDF we obtain 
as{Mz) = 0.127 which is about 15% larger than the value of 0.110 in CTEQ2M PDF 
P7| . We found that if we rescale the ti production rates obtained from CTEQ2L 
PDF with born level amplitudes by the ratio of a'^{Q, Aqcd) from CTEQ2M and 
that from CTEQ2L, which yields 0.7 for Q = Mz, then our total rates are in good 
agreement with those obtained using NLO PDF and NLO amplitudes see, for 



example. Reference [^. Hereafter we shall use the scaled results for our rates. The 
constituent cross sections are all calculated at tree level for simplicity to study the 
kinematics of the top quark and its decay products. 

To include the production rates for both single-t and single-t events at pp colliders. 



^ In reality, the integrated luminosity can be higher than the ones used here. For instance, with 
a couple of years of running a 2 TeV Tevatron can accumulate, say, 10fb~^ luminosity. Similarly, it 
is not out of question to have a 4TeV Di-TeV to deliver an integrated luminosity of about 100 fb~^. 
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Figure 3.2: Rate in [pb] for qq, gg — > tt, q'g{W^g) —>■ qtb, q'q —>■ W* —>■ tb and 
gb — > W~t at various energies of pp colliders. 
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Figure 3.3: Rate in [pb] for qq, gg tt, q'g{W^g) 
gb — > W~t at various energies of pp colliders. 
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Table 3.1: Rates of the above processes for rrit = 180(140) GeV. (Branching ratios 
are not included here.) For y/S = 2 TeV and 4 TeV we include rates for a pp machine. 
At = 14 TeV the rates are for a pp machine. For the single-top rates we only 
include single-t production. 





Cross Section (pb) 




QQ, 99 tt 


q'g — > qtb (or q'b qt) 


q'q ^W* ^tb 


gb ^ W-t 


2 


4.5(16) 


1(2) 


0.3(0.8) 


0.1(0.3) 


4 


26(88) 


7(11) 


0.8(2.1) 


1.3(2.9) 


14 


430(1300) 


100(140) 


4.6(11) 


3.6(8.8) 



a factor of 2 should be multiplied to the single-t rates shown in Figures p.2| and p.3| 
because the parton luminosity for single-t production is the same as that for single-t. 
Similarly, at pp colliders the rates should be multiplied by about 1.5 for the center- 

of-mass energy (v^) of the collider up to approximately 4 TeV, but almost a factor 

of two at higher energies (say, > 8 TeV up to about 14 TeV) because the relevant 
parton luminosities for producing a single-t and a single-t event in pp collisions are 



different. As shown in Figures and |3.3| the total rate for single-top production is 

about the same at pp and pp colliders for > 8 TeV because the relevant valence 
and sea quark parton distributions are about equal for 100 GeV < rrit < 300 GeV. 

For smaller y/S, up to about 4 TeV, a pp collider is preferred over a pp collider for 
heavy top quark production because of its larger parton luminosities. Similarly, for ti 



pair productions at small y S, the quark initiated process qq tt is more important 

than the gluon fusion process gg — > ti. At y/S from 8 to 14 TeV the ti rate is about 
the same in pp and pp collisions because the gg ti subprocess becomes dominant. 

For later reference in this paper, we show the rates of the above processes in 
Table ^]1] for rrit = 180(140) GeV. (Branching ratios are not included here.) For 



y/S = 2 and 4 TeV we include only the rates for a pp machine, whereas at y/S = 14 
TeV the rates are for a pp machine. Again, for the single-top rates we only include t 
production. 

Both in Figures and |3.3| and Table |3.1| , we have given the cross section of 
single-top quarks produced from either the qg{W^g) — > qtb or q'b qt processes. 
From now on, we will refer to this production rate as the rate of the VT-gluon fusion 
process. The single-top quark produced from the H^-gluon fusion process involves a 
very important and not yet well-developed technique of handling the kinematics of 
a heavy b parton inside a hadron. Thus the kinematics of the top quark produced 
from this process can not be accurately calculated yet. However, the total event rate 
for single-top quark production via this process can be estimated using the method 
proposed in Reference The total rate for the VT-gluon fusion process involves the 
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Figure 3.4: Feynman diagrams illustrating the subtraction procedure for calculating 
the total rate for l^-gluon fusion: q'b ^ qt ® q'g(W~^g) — ^ qtb Q {g ^ bb® q'b qt). 

0{a^) (2 2) process q'b — > qt plus the 0{a^as) (2 3) process q'giyV'^g) qtb 
(where the gluon splits to bb) minus the splitting piece g ^ bb ® q'b qt in which bb 
are nearly coUinear. These processes are shown diagrammatically in Figure The 
helicity amplitudes and the cross sections for these processes are given in Appendices 
A and B respectively. 

The splitting piece is subtracted to avoid double counting the regime in which 
the b propagator in the (2 3) process is close to on-shell.0 The procedure is to 
resum the large logarithm Q;sln(m^/m^) in the VT-gluon fusion process to all orders 
in as and include part of the higher order Oi^oP'as) corrections to its production rate. 
{rrih is the mass of the bottom quark.) We note that to obtain the complete Oi^a^as) 
corrections beyond just the leading log contributions one should also include virtual 
corrections to the (2 — 2) process, but we shall ignore these non-leading contributions 
in this work. Using the prescription described as above we found that the total rate 
of the l^-gluon fusion process is about 25% less as compared to the (2 2) event 
rate for = 180 (140) GeV regardless of the energy or the type {i.e., pp or pp) of 



the machine. In Figures ^]5| and |3.6| we show the total rate of ly-gluon fusion versus 
rrit with scale Q = mj as well as a breakdown of the contributing processes at the 
Tevatron, the Di-TeV and the LHC. 



To estimate the uncertainty in the production rate due to the choice of the scale 
Q in evaluating the strong coupling constant as and the parton distributions, we 



show in Figure the scale dependence of the W^-gluon fusion rate. As shown in 
the figure, although the individual rate from either (2 — 2), (2 ^ 3) or the splitting 
piece is relatively sensitive to the choice of the scale, the total rate as defined by 
(2 ^ 2) + (2 ^ 3) - (splitting piece) only varies by about 30% for Mvk/2 <Q < 2mt 
at the Tevatron. At the Di-TeV and the LHC, it varies by about 30% and 10%, 
respectively. Based upon the results shown in Figure |^, we argue that Q < Mw/2 
probably is not a good choice as the relevant scale for the production of the top 
quark from the VT-gluon fusion process because the total rate rapidly increases by 
about a factor of 2 in the low Q regime. In view of the prescription adopted in 

^ The total rate of the (2 3) process is extensively discussed in the Appendix C. 
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Figure 3.5: Rate in [pb] for single-t production: q'b — > (2 — > 2), q'g — > qtb (2 — > 3) 
and the splitting piece g ^ bb ® q'b ^ qt in which bb are coUinear. The rates are for 
pp coUiders. 
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Figure 3.6: Rate in [pb] for single-t production: q'b — > (2 — > 2), q'g — > qtb (2 — : 
and the splitting piece g ^ bb ® q'b ^ qt in which bb are coUinear. The rates are 
pp coUiders. 
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Figure 3.7: Rate of PF~gluon fusion process versus scale Q for mt = 180 GeV and 
v^ = 2TeV. 
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calculating the total rate, the only relevant scales are the top quark mass rrit and 
the virtuality of the ly-line in the scattering amplitudes. Since the typical transverse 
momentum of the quark (g), which comes from the initial quark {q') after emitting 
the VT-line, is about half of the l^-boson mass, the typical virtuality of the In- 
line is about Mvi//2 ^ 40GeV. rrib — 5 GeV is thus not an appropriate scale to be 
used in calculating the W^-gluon fusion rate using our prescription. We note that 
in the (2 — 2) process the b quark distribution effectively contains sums to order 
[a^ ln((5/'nib)]" from n-fold coUinear gluon emission, whereas the subtraction term 
(namely, the splitting piece) contains only first order in asln{Q /mi,). Therefore, as 
Q ^ nih the (2 2) process picks up only the leading order in asln{Q/mi,) and 
so gets largely cancelled in calculating the total rate. Consequently, as shown in 
Figure ^.7[ the total rate is about the same as the (2 — > 3) rate for Q rrif,. We 
also note that at Q — Mwf^, the (2 — > 2) and (2 3) processes have about the 
same rate. As Q increases the (2 — ^ 2) rate gradually increases while the (2 — > 3) 
rate decreases such that the total rate is not sensitive to the scale Q. It is easy to see 
also that the total rates calculated via this prescription will not be sensitive to the 
choice of PDF although each individual piece can have different results from different 
PDF's, based upon the factorization of the QCD theory ^0 . 



Another single-top quark production mechanism is the Drell-Yan type process 
q'q W* —>■ th. As shown in Figures |3.2| and p.3| , for top quarks with mass on the 
order of 180 GeV the rate for W* production is about one fifth that of ly-gluon fusion 

at ^fS = 2^eV. The W* process becomes much less important for a heavier top quark. 
This js because at higher invariant masses s (for producing a heavier top quark) of 
the tb system, W* production suffers the usual 1/s suppression in the constituent 
cross section. However, in the ly-gluon fusion process the constituent cross section 
does not fall off as 1/s but flattens out asymptotically to 1/M^. (The analytical 
results of these amplitudes are given in Appendix A for reference.) For colliders 
with higher energies, therefore with large range of s, the W* production mechanism 
for heavy top quarks becomes much less important. However, the kinematics of 
the top quarks produced from this process are different from those in the W^-gluon 
fusion events. Moreover, possible new physics may introduce a high mass state (say, 
particle V) to couple strongly with the tb system such that the production rate from 
q'q W* —>■ V —>■ tb can largely deviate from the SM W* rate.[| We will however 
not discuss it in detail here because its rate is highly model dependent. 

The ly-gluon fusion process becomes more important for a heavier top quark. 
Why? Effectively, the M^-gluon fusion process can be viewed as the scattering of a 
longitudinal l^-boson (Wl) with gluon to produce a top quark and a bottom anti- 
quark (W^g tb) after applying the effective- approximation [^. For large s 
this scattering process is equivalent to {4>'^g — *■ tb) where 0"'" is the corresponding 
Goldstone boson of the gauge boson due to the Goldstone Equivalence Theorem 



|33| , 0. Since the coupling of t-b-(j) is proportional to the mass of the top quark, 
the constituent cross section of the IV-gluon fusion process grows like m^/M^ when 
rrit increases. This explains why the H^-gluon fusion rate only decreases slightly as 
the mass of the top quark increases even though both the parton luminosity and the 
available phase space decrease for a heavier top quark. In contrast, the tt pair produc- 
tion rate from the QCD processes decreases more rapidly as rrit increases because the 



^ This is similar to the speculations made in Reference ||3l|| for having some high mass resonance 
in tt production. 
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constituent cross section of qq, gg ti goes as 1/s and the phase space for producing 
a tt pair is smaller than that for producing a single-t. Therefore, the H^-gluon fusion 
process becomes more important for the production of a heavy top quark. 

Before closing this chapter, we note that the Effective-iy approximation has been 
the essential tool used in studying the strongly interacting longitudinal W system 
to probe the symmetry breaking sector at the supercolliders such as the LHC |p5|| . 
By studying single-top production from the W^-gluon fusion process at the Teva- 
tron, one can learn about the validity of the Effective-H^ approximation prior to the 
supercolliders. 



Chapter 4 

Measuring the Top Quark Mass 



By the year 2000, we expect results from the Tevatron (with Ifb^^) and results 
from LEP-200, giving an error of about 50MeV on Mw- Due to Veltman's screening 
theorem, the low energy data are not sensitive to the mass of the Higgs boson . For 
a heavy Higgs boson, the low energy data can at most depend on ruH logarithmically 
up to the one loop level. Therefore, within the SM one needs to also know the mass 
of the top quark to within about 5 GeV to start getting useful information on rriH 
with an uncertainty less than a few hundred GeV. This can be done by studying 
radiative corrections to the low energy data which include LEP, SLC, and neutrino 

m 



experiments 
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(Of course, rriH will be measured to better precision if 



it is detected from direct production at colliders.) 

How accurate can the mass of the top quark be measured at hadron colliders? At 
hadron colliders, rrit can be measured in the tt events by several methods [22, 37|. 
The first method is to use the lepton+jet decay mode of the ti pair, as shown in 



Figure |4.1| . This is done by reconstructing the invariant mass of the three jets in 
the opposite hemisphere from the isolated lepton i (= eor^u) in t — > bW{—> iu), 
and requiring that two of the three jets reconstruct to a. W and the third be tagged 
as a 6-jet. The second method is to use the di-lepton decay mode of the ti pair. 



as shown in Figure O. Here it is required that both W^s decay leptonically. In 
addition, one of the 6's must decay semileptonically to measure the mass distribution 
of the non-isolated lepton ib (from b decay) and one of the two isolated leptons (£i 
and £2 from decay) which is closer to ib- The third method is to measure the 
cross section of the di-lepton decay mode of the ti pair. At the LHC, there will be 
about 10^ ti pairs produced in one year of running for nit < 200 GeV. With such a 
large number of events, the ATLAS and CMS collaborations concluded that can 
be measured with a precision of < 5 GeV using the first method described above and 
with about a factor of 2 improvement using the second method pSl E^. A similar 



conclusion was also drawn by the CDF and the D0 collaborations for the Tevatron 
with Main Injector after the upgrade of their detectors [Q. This is remarkable given 
that the ti cross section at the Tevatron is smaller by about two orders of magnitude 
as compared with that at the LHC, as shown in Figures p.2| and p.3|. 



Next, we would like to discuss how to measure the mass of the top quark in 
the VT-gluon fusion process. Since rrit has been measured by the FNAL CDF and 
D0 groups in the ti events p3| , p!6[| , why do we care? To check whether it is a SM top 
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quark, we should verify its production rate predicted by the SM for other production 
processes such as the single-top quark process. 

Suppose the coupling of t-b-W is not of the SM nature, then we would find that 
the single-top quark production rate of the VT-gluon fusion process is different from 
the SM prediction because its production rate is directly proportional to the square 
of this coupling. (We will discuss more on this point in Chapter 6.) Hence, without 
knowing the nature of the t-b-W interactions one can not use the production rates of 
the single-top quark events to measure m^. Alternatively, we propose two methods to 
measure rrit in the single-top quark events. We will refer to them as the fourth and the 
fifth method. The fourth method is a slight variation of the second method. Instead 
of measuring the invariant mass of the leptons, we propose to directly measure the 
invariant mass {rribe) of the £ and b in t bW{—^ iu). We expect that the efficiency 
of b tagging using the displaced vertex is higher for detecting a heavier top quark, 
and the b jet energy measurement is better for b having larger transverse momentum 
from a heavy top quark decay. Thus mbe can be used to measure the mass of a 
SM top quark. The details of our Monte Carlo study are given in Chapter 8 for a 
single-top quark event. In the tt event there are two 6's, therefore this method may 
not work as well as in the single-top event which only contains one b. However it is 
not entirely impossible to use this method because, as shown in Figure |4.3| , the sum 
of the invariant mass distributions of bi and bi for a 180 GeV top quark still show 
a bump near the region that the distribution of rribe peaks. (With a larger sample 
of ti events one might be able to afford using the electric charge of the soft-lepton 
from 6-decay to separate b from b on an event-by-event basis at the cost of the small 
branching ratio of 6 ^ /i + X, of about 10%.) We will explain in more detail how 
to use /Long (the fraction of longitudinal H^-bosons from top quark decay), derived 
from the distribution of mf,e, to measure in Chapter 6. 

The fifth method is to reconstruct the invariant mass of the top quark in the 
t —>■ bW{^ iu) decay mode by measuring the missing transverse momentum and 
choosing a two-fold solution of the longitudinal momentum of the neutrino from the 
mass constraint of the W boson. In Chapter 8 we conclude that it is possible to 
measure rrit using either of these last two methods to a precision of 5 GeV at the 

Tevatron (\^ = 2TeV) with lfb~^ integrated luminosity. We also find that after 
applying all the kinematical cuts to suppress the dominant background W + bb, at 
most 10% of W* events contribute to single-top production for a 180 (140) GeV top 
quark. The SM W* production rate is already much smaller than the 14^-gluon fusion 
rate for a heavier top quark, therefore the contribution from the W* is not important 
in our study although we do include its small effects in our analysis as described in 
Chapter 8. 
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Figure 4.3: Distributions of rribi (solid) and mg^ (dash) in tt events for a 180 GeV 
top quark. 



Chapter 5 

Measuring the Top Quark Width 



As shown in Reference |^T| the intrinsic width of the top quark can not be measured 
at a high energy hadron colhder such as the LHC through the usual QCD processes.^] 
For instance, the intrinsic width of a 150 GeV Standard Model top quark is about 1 
GeV, and the full width at half maximum of the reconstructed top quark invariant 
mass (from t — >• bW{^ jets) decay mode) is about 10 GeV after including the 
detector resolution effects by smearing the final state parton momenta. Here, the 
ratio of the measured width and the intrinsic width for a 150 GeV top quark is about 
a factor of 10. For a heavier top quark, this ratio may be slightly improved because 
the jet energy can be better measured. (The detector resolution AE/E for a QCD jet 



with energy E is proportional to l/yE.) A similar conclusion was also given from a 
hadron level analysis presented in the SDC Technical Design Report which concluded 
that reconstructing the top quark invariant mass gave a width of 9 GeV for a 150 
GeV top quark Is there a way to measure the top quark width T(t ^ bW^ 



say, within a factor of 2 or better, at hadron colliders? Yes, it can in principle be 
measured in single-top events. 

The width r(t — > bW~^) can be measured by counting the production rate of top 
quarks from the W-b fusion process which is equivalent to the IV-gluon fusion process 
by a proper treatment of the bottom quark and the W boson as partons inside the 
hadron. The W^-boson which interacts with the 6-quark to produce the top quark can 



be treated as an on-shell boson in the leading log approximation The result 

is that even under the approximations considered, a factor of 2 uncertainty in the 
production rate for this process gives a factor of 2 uncertainty in the measurement 
of r(i(: — > bW^). This is already much better than what can be measured from the 
invariant mass distribution of the jets from the decay of top quarks in the tt events 
produced via the usual QCD processes. More precisely, as argued in Chapter 3, the 
production rate of single-top events at the Tevatron can probably be known within 
about 30%, thus it implies T(t — ^ bW~^) can be measured to about the same accuracy.0 

^ In Reference Q], the effects of QCD radiation in top quark decay (at one loop level) to the 
measurement of toj in tt events produced in hadron collisions was studied. It was concluded that 
the peak position of the mt distribution remains about the same as the tree level result, but the 
shape is different. It was also found that the mi,e distribution is not sensitive to QCD radiation in 
top decay. 

^ Strictly speaking, from the production rate of single-top events, one measures the sum of all the 
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Therefore, this is an extremely important measurement because it directly tests the 
couplings of t-b-W . 

H^-gluon fusion can also tell us about the CKM matrix element \Vtb\- Assuming 
only three generations of quarks, the constraints from low energy data together with 
unitarity of the CKM matrix require \Vtb\ to be in 0.9988 to 0.9995 at the 90% 
confidence level 0. As noted in Reference P the low energy data do not preclude 
there being more than three generations of quarks (assuming the same interactions 
as described by the SM). Moreover, the entries deduced from unitarity might be 
altered when the CKM matrix is expanded to accommodate more generations. When 
there are more than three generations the allowed ranges (at 90% CL) of the matrix 
element \Vtb\ can be anywhere between and 0.9995 Since \Vtb\ is directly involved 
in single-top production via H^-gluon fusion, any deviation from SM value in \ Vtb\ will 
either enhance or suppress the production rate of single-top events. It can therefore be 
measured by simply counting the single-top event rates. For instance, if the single-top 
production rate is measured to within 30%, then \Vtb\ is determined to within 15%. 

In conclusion, after the top quark is found, the branching ratio of t — >■ bW~^{^ 
C.'^v) can be measured from the ratio of (2£+ jets) and {11+ jets) rates in ti events. 
The measured single-top quark event rate is equal to the single-top production rate 
multiplied by the branching ratio of t — > bW~^{—>- l^v) for the (1£ + jets) mode 
and the same t-b-W couplings appearing in the decay of t in this process appear 
also in the production of t. Thus, a model independent measurement of the decay 
width r(t — > bW~^) can be made by simply counting the production rate of t in the 
IV-gluon fusion process. Should the top quark width be found to be different from 
the SM expectations, we would then have to look for non-standard decay modes of 
the top quark. We note that it is important to measure at least one partial width 
(say, r{t —>■ bW~^)) precisely in order to discriminate between different models of new 
physics, if any. In the SM, the partial width r{t —>■ bW~^) is about the same as the 
total width of the top quark at the tree level because of the small CKM matrix element 
\Vts\, thus measuring the single-top quark production rate measures the lifetime of 
the top quark. 



possible partial decay widths, such as r{t bW^) + T{t sW^) + T{t dW^) + • ■ •, therefore, 
this measureraent is really measuring the width of T(t XW^) where X can be more than one 
particle state as long as it originates from the partons inside the proton (or anti-proton). In the SM, 
r{t ^ bW'^) is about equal to the total width of the top quark. 



Chapter 6 



Top Quark Couplings to the W 
Gauge Boson 



It is equally important to ask what kind of interactions the t-b-W vertex might involve 
4^ . For instance, one should examine the form factors of t-b-W which result from 
higher order corrections due to SM strong and/or electroweak interactions. It is 
even more interesting to examine these form factors to test the plausibility of having 
nonuniversal gauge couplings of t-b-W due to some dynamical symmetry breaking 
scenario |19|. 

The QCD and the electroweak |^ corrections to the decay process t — > bW^ 
in the SM have been done in the literature. The most general operators for this 
coupling are described by the interaction lagrangian 



L = ^ 



9_ 
V2 



1 



■d.W-ba^''{f^P^ + /«P+)t 



w 



9 

'V2 



M, 



w 



(6.1^ 



where P± 



^1 ±75), ia^ 



'ib'^'T^'] superscript * denotes the complex 



conjugate. In general, the form factors fi'^ and /;f can be complex. Note that in 
Equation (|6.1|), if there is a relative phase between ff" and or between ff^ and 
then CP is violated. For instance, in the limit of = , a CP-violating observable 
will have a coefficient proportional to lm{fif^ ) for a left-handed bottom quark. 



and Im(//^/^ ) for a right-handed bottom quark (We will discuss more on CP 
violation in Chapter 7.) If the VT-boson can be off-shell then there are additional 
form factors such as 



d'-wM^p. + fip+)t + d^wnuf P- + p+)b 



(6.2) 



which vanish for an on-shell W^-boson or when the off-shell VT-boson couples to 
massless on-shell fermions. Here, we only consider on-shell H^-bosons for mt > 
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Mw + rrif,. At tree level in the SM the form factors are = 1 and = f2 = fi^ = ^■ 
These form factors will in general affect the experimental observables related to the 
top quark, such as the fraction of longitudinal W's produced in top quark decays. 

The fraction (/Long) of longitudinally polarized ly-bosons, produced in the rest 

frame of the decaying top quark, strongly depends on the form factors fi'^ and 

/^'^, as shown in Appendix C. Hence, /Long is a useful observable for measuring 
these form factors. The definition of /Long is simply the ratio of the number of 
longitudinally polarized H^-bosons produced with respect to the total number of W- 
bosons produced in top quark decays: 

f n^w = 0) 

r(Aw = 0) + = -)+ nxw = +) • ^ ^ ^ 

We use T{\]y) to refer to the decay rate for a top quark to decay into a ly-boson with 
polarization A^y. (Av^^ = — , +, denotes a left-handed, right-handed, and longitudinal 
ly-boson.) Clearly, the polarization of the ly-boson depends on the form factors /i 
and /2.Q Therefore, one can measure the polarization of the ly-boson to measure 
these form factors. As shown in Appendix C, the polarization of the W^-boson can be 
determined by the angular distribution of the lepton, say, e"*" in the rest frame of in 
the decay mode t — >■ bW~^{^ e^v). However, the reconstruction of the W-hoson rest 
frame (to measure its polarization) could be a non-trivial matter due to the missing 
longitudinal momentum {Pz) (with a two-fold ambiguity) of the neutrino (z/) from W 



decay. Fortunately, as shown in Equation ( |6.4| ), one can determine the polarization 
of the W-hoson without reconstructing its rest frame by using the Lorentz-invariant 
observable mbe, the invariant mass of h and e from t decay. 

The polar angle 6'*+ distribution of the e"*" in the rest frame of the boson, 
whose z-axis is defined to be the moving direction of the boson in the rest frame 
of the top quark, can be written in terms of m^e through the following derivation: 

EeEb -pe-Pb 



cos 



iPellPfel 

Pe-Pb ^ ^ _ 2mgg 



(6.4) 



The energies E^, and Ei, are evaluated in the rest frame of the boson from the 
top quark decay and are given by 

Ee - jrrj , \X)e\ - \l Et - mt, 

MVlwr 



ml - - ml 



^ The fraction of longitudinal Ws in top quark decays contributed from the form factor //^ is 
the same as that from f[ |44[. 
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me (m^) denotes the mass of e"*" (z/g) for the sake of bookkeeping. The first hne 
in Equation ( |6.4| ) is exact when using Equation (|6.5|) , while the second hne of Equa- 
tion ( |6.4| ) holds in the limit of rrih = 0. It is now trivial to find /Long by first calculating 
the cos distribution then fitting it according to the decay amplitudes of the W- 
boson from top quark decay, as given in the Appendix C. In what follows we will 
show how to use the distribution of m^e to measure the mass of the top quark and its 
couplings to the VT-boson. 

In Reference [0, we considered the effective couplings 

1 CC 

W-t,-b,: ^^^7m(1 - 75) (6.6) 

and 

cc 

W-t^-bn: ^^7.(1 + 75) (6.7) 

derived from an electroweak chiral lagrangian with the symmetry SU(2)^ x U(l)y 

broken down to U(1)em- (Here, = /f' - 1, and k%'^ = f^.) At the Tevatron 
and the LHC, heavy top quarks are predominantly produced from the QCD process 
gg, qq ti and the VT-gluon fusion process qg{Wg) — > tb, ib. In the former process, 
one can probe and kr'"'" from the decay of the top quark to a bottom quark and 
a VT-boson. In the latter process, these non-standard couplings can also be measured 
by simply counting the production rates of signal events with a single t or t. Let us 
discuss kl*"'" and k^j'-^'" in more detail as follows. 



6.1 From the Decay of Top Quarks 

To probe and kj?*"'" from the decay of the top quark to a bottom quark and a TV- 
boson, one needs to measure the polarization of the IV-boson, which can be measured 
from the distribution of the invariant mass rribe. For a massless b, the IV-boson from 
top quark decay can only be either longitudinally or left-hand polarized for a left- 
hand charged current {kr"^ = 0). For a right-hand charged current {kl"'" = — 1) 
the IV-boson can only be either longitudinally or right-hand polarized. (Note that 
the handedness of the Vl^-boson is reversed for a massless b from i decays.) This 
is a consequence of helicity conservation, as diagrammatically shown in Figures |6.1| 



and |6.2| for a polarized top quark. In these figures we show the preferred moving 
direction of the lepton from a polarized W^-boson in the rest frame of a polarized 
top quark for either a left-handed or a right-handed t-b-W vertex. As indicated in 
these figures, the invariant mass mM depends on the polarization of the IV-boson 
from the decay of a polarized top quark. Also, mM is preferentially larger for a 
pure right-handed t-b-W vertex than a pure left-handed one. This is clearly shown 



in Figure |6.3| , in which the peak of the tum distribution is shifted to the right and the 
distribution falls off sharply at the upper mass limit for a pure right-handed t-b-W 
vertex. In terms of cos^|, their difference is shown in Figure |6.4|. However, in both 



cases the fraction (/Long) of longitudinal W^s from top quark decay is enhanced by 
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Figure 6.1: For a left-handed t-b-W vertex. 
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Figure 6.2: For a right-handed t-h-W vertex. 
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Figure 6.3: ruhi distribution for SM top quark (solid) and for pure right- 
handed t-b-W couphng of t6Vr(dash). 
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Figure 6.4: cos 9} distribution for SM top quark (solid) and for pure right- 
hand t-b-W couphng of i6VF(dash). 
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rrit^ /2My/'^ as compared to the fraction of transversely polarized VT's [P], namely, 



2M''' 

/Long = - ;;;2~ • (6. 
1 + 



2M,i 



Therefore, for a heavier top quark, it is more difficult to untangle the and k,r"^ 
contributions. On the other hand, because of the very same reason, the mass of a 
heavy top quark can be accurately measured from /Long (discussed below) irrespective 
of the nature of the t-b-W couplings (either left-handed or right-handed). 

The QCD production rate of ti is obviously independent of kl^'" and kr"'" . (Here 
we assume the electroweak production rate of qq A, Z ti remains small as in 
the SM.) Let us estimate how well the couplings kl^'" and kr"^ can be measured at 
the Tevatron, the Di-TeV, and the LHC First, we need to know the production rates 
of the top quark pairs from the QCD processes. As shown in Table |3.1| , the QCD 
production rate of gg,qq ti for a 180 GeV top quark is about 4.5 pb, 26 pb and 
430 pb at the Tevatron, the Di-TeV, and the LHC, respectively. For simplicity, let's 
consider the £^ + > 3 jet decay mode whose branching ratio is Br = 2|| = ^, where 

the charged lepton can be either or n^. We assume the experimental detection 
efficiency (e), which includes both the kinematic acceptance and the efficiency of b- 
tagging, to be 15% for the signal event Let's further assume that there is no 

ambiguity in picking up the right b (b) to combine with the charged lepton 

to reconstruct t (or i), then in total there are 4.5 pb x lO'^pb"^ x ^ x 0.15 = 200 

reconstructed ti events to be used in measuring kl and k,r at VS = 2 TeV. The 
same calculation at the Di-TeV and the LHC yields 1100 and 19000 reconstructed tt 
events, respectively. Given the number of reconstructed top quark events, one can fit 
the tum distribution to measure kl'"'" and hr"'~^ . For example we have done a study 
for the Tevatron. Let us assume the effects of new physics only modify the SM results 

[f^ = 1 and = at Born level) slightly and the form factors /^'^ are as small 
as expected from the usual dimensional analysis [^, |5^ .| We summarize our results 



on the accuracy of measuring fi'^ for various luminosities in Table |6.1| ||51|| . (Only 
statistical errors are included at the 95% confidence level.) 



In the same table (i.e., Table |6.1| ) we also show our estimate on how well the mass 
of the top quark can be measured from /Long- By definition of /Long, for a SM top 

quark {i.e., fi=l and ff = 0), the distribution of cos 6*1 has the functional form 



F(coso;) ~ (1 - (i^) + A„„ (!^' 



(6,9) 



9 L R 

The coefficients of the form factors /2 ' , assumed to be induced through loop effects, wiU be 
a factor of -r-^ smaller than that of the form factors A ' . 
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Table 6.1: Results on the accuracy of measuring //"'^ for various luminosities. (Only 
statistical errors are included at the 95% confidence level.) 



Integrated 


Number of 








Luminosity 

fb-i 


reconstructed 
tt events. 


ft 




Amt 

mt 


1 


200 


8% 


±0.5 


4% 


3 


600 


4% 


±0.3 


2% 


10 


2000 


2% 


±0.2 


1% 



Therefore, /Long can be calculated by fitting with the distribution of cos6'|, or equiv- 

alently with the distribution of rribe- We prefer to measure kl'"'" and k_r'"'" using the 
distributions of ttim than of cos^^ because the former can be directly calculated from 
the measured momenta of b and i. However, to convert from the distributions of rribe 
to cosO^, as given in Equation ( |6.4|) , the effects from the width of the M^-boson and 
the top quark might slightly distort the distribution of cos^^. (Notice that in the 
full calculation of the scattering amplitudes the widths of the ly-boson and the top 
quark have to be included in the Breit-Wigner form to generate a finite event rate.) 

However, in reality, the momenta of the bottom quark and the charged lepton will 
be smeared by detector effects and another problem in this analysis is the identifica- 
tion of the right b to reconstruct t. There are three possible strategies to improve the 
efficiency of identifying the right b. One is to demand a large invariant mass of the tt 
system so that t is boosted and its decay products are collimated. Namely, the right 
b will be moving closer to the lepton from t decay. This can be easily enforced by 
demanding leptons with a larger transverse momentum. Another is to identify the 
soft (non-isolated) lepton from b decay (with a branching ratio Br (6 — > fJ^^X) ~ 10%). 
The other is to statistically determine the electric charge of the &-jet (or 6-jet) to be 
1/3 (or —1/3) All of these methods may further reduce the reconstructed signal 
rate by an order of magnitude. How will these affect our conclusion on the determi- 
nation of the non-universal couplings kl'"'" and i^r"'"'! It can only be answered by 
detailed Monte Carlo studies which are yet to be done. 



6.2 From the Production of Top Quarks 

Here we propose another method to measure the couplings kl*"'" and kr"'" from the 
production rate of the single-top quark process. 

For rrit = 180 GeV, the sum of the production rates of single-t and single-t events 

is about 2 pb and 14 pb for = 2 TeV and = 4 TeV respectively. The branching 

ratio of interest is Br = |. The kinematic acceptance of this event at = 2 TeV 
is about 0.55, as shown in Chapter 8. Assuming the efficiency of 6-tagging is about 
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30%, then there will be 2 pb x 10^ pb ^ x | x 0.55 x 0.3 = 75 events reconstructed 

for a Ifb^^ integrated luminosity. At = 4TeV, as shown in Chapter 8, the 
kinematic acceptance of this event is about 0.40 which, from the above calculation, 
yields about 3700 reconstructed events for 10fb~^ integrated luminosity. Based on 
statistical error alone, this corresponds to a 12% and 2% measurement on the single- 
top cross section. A factor of 10 increase in the luminosity of the collider can improve 
the measurement by a factor of 3 statistically. Taking into account the theoretical 
uncertainties, as discussed in Chapter 3, we examine two scenarios: 20% and 50% 
error on the measurement of the cross section for single-top production. The results, 
which are not sensitive to the energies of the colliders considered here (either 2 TeV 
or 4 TeV), are shown in Figure for a 180 GeV top quark at the Tevatron. We 
found that ki!"^ and n^^^ are well constrained inside the region bounded by two 
(approximate) ellipses (cf. Appendix A). To further determine the sizes of kl"'" and 
i^R^'~^ one needs to study the kinematics of the decay products, such as the charged 
lepton £, of the top quark. Since the top quark produced from the VT-gluon fusion 
process is almost one hundred percent left-hand (right-hand) polarized for a left- 
hand (right-hand) t-b-W vertex, the charged lepton from t decay has a harder 
momentum for a right-handed t-b-W coupling than for a left-handed coupling. (Note 
that the couplings of light-fermions to VT-boson have been well tested from the low 
energy data to be left-handed as described in the SM.) As shown in Figures ^]l] 
and |6.2| , this difference becomes smaller when the top quark is much heavier because 
the l^-boson from the top quark decay tends to be more longitudinally polarized. 

A right-hand charged current is absent in a linearly SU(2)j^ invariant gauge theory 

with massless bottom quark. In this case, kr"'" = 0, then kl'"'" can be constrained 
to within about -0.08 < k^^^ < 0.03 (-0.20 < kl^^ < 0.08) with a 20% (50%) 
measurement on the production rate of single-top quark at the Tevatron [jl9l . (Here 



we assume the experimental data agrees with the SM prediction within 20% (50%).) 
This means that if we interpret (1 -|- i^lP'") as the CKM matrix element |Vjfe|, then 
|Vtf,| can be bounded as \Vtb\ > 0.9 (or 0.75) for a 20% (or 50%) measurement on the 
single-top production rate. 

Before closing this chapter, we remark that in the Refs. |jl9[ and [|53[ some bounds 
on the couplings of kl'"'" and k^j'-^'" were obtained by studying the low energy data 
with the assumption that the effects of new physics at low energy can only modify 
the couplings of kl'^'^ and kr'"'^ but not introduce any other light fields in the ef- 
fective theory. However, nature might not behave exactly in this way. It is possible 
that some light fields may exist just below the TeV scale, then the bounds obtained 
from Refs. |TPI] and |E3I] may no longer hold. Thus, it is important to have direct 



measurements on all the form factors listed in Equation ( |6.1[ ) from the production of 
top quarks, in spite of the present bounds on k's derived from radiative corrections 
to low energy data. 
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Figure 6.5: Constraint on \kl~'^\ and kr-"^ given 20% and 50% error in measurement 
of Standard Model rate for 1^-gluon fusion. Curves are identical for mt — 140 GeV 
and mt = 180 GeV. 



Chapter 7 

Probing CP Properties in Top 
Quarks 



It is known that explicit CP violation requires the presence of both the CP non- 
conserving vertex and the complex structure of the physical amplitude. Due to the 
origin of this complex structure, the possible CP-violating observables can be sepa- 
rated into two categories. In the first category, this complex structure comes from 
the absorptive part of the amplitude due to the final state interactions. In the sec- 
ond category, this complex structure does not arise from the absorptive phase but 
from the correlations in the kinematics of the initial and final state particles involved 
in the physical process. Hence, it must involve a triple product correlation [i.e., a 
Levi-Civita tensor). 

To distinguish the symmetry properties between these two cases, we introduce 
the transformation T, as defined in Reference [|^, which is simply the application of 



time reversal to all momenta and spins without interchanging initial and final states. 

The CP-violating observables in the first category are CP-odd and CPT-odd, while 

those in the second category are CP-odd and CPT-even. Of course, both of them are 
CPT-even. 

To illustrate the above two categories, we consider CP-violating observables for 
the decay of the top quark. Consider the partial rate asymmetry 

^ r(t-^bw+)-Tit^bw-) 

-^'^ - T{t^bW+) + T{i-^bW-y ^ ' 

A non- vanishing Abw clearly violates CP and CPT, therefore this observable belongs 
to the first category. We note that because of CPT invariance, the total decay width 
of the top quark r(t) has to equal the total decay width of the top anti-quark T(t). 
Thus, any non-zero Abw implies that there exists a state (or perhaps more than one 
state) X such that t can decay into X and t into X. The absorptive phase of t — > bW~^ 
is therefore generated by re-scattering through state X, i.e., t ^ X ^ bW^, where 
X 7^ bW~^ because the final state interaction should be off-diagonal . 
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Next, let's consider the observable of the second category. In the decay of t ^ 
bW~^{^ i'^i'e), for a polarized t quark, time reversal invariance (T) is violated if the 
expectation value of 

^ty< Pb- Pe+ (7.2) 



is not zero Assuming CPT invariance, this implies CP is violated. Therefore, this 



observable is CP-odd but CPT-even. A non- vanishing triple product observable, such 
as that in Equation (|7.2|), from the decay of the top quark violates T. However, it may 
be entirely due to final state interaction effects without involving any CP-violating 
vertex. To construct a truly CP-violating observable, one must combine information 
from both the t and t quarks. For instance, the difference in the expectation values 
of eft X Pf Pi+ and (f^ x p^- p^- would be a true measure of an intrinsic CP violation. 

There have been many studies on how to measure the CP-violating effects in the 
ti system produced in either electron or hadron collisions. (For a review, see a recent 
paper in Reference ||56|.) At hadron colliders, the number of ti events needed to 



measure a CP-violating effect of the order of 10~^ — 10~^ is about 10^ — 10^. To 
examine the potential of various current and future hadron colliders in measuring 
the CP-violating asymmetries, we estimate the total event rates of ti pairs for a 180 
GeV SM top quark produced at these colliders. At the Tevatron, the Di-TeV and the 
LHC, an integrated luminosity of 10, 100 and 100 fb~^ will produce about 4.5 x 10^, 
2.6 X 10^ and 4.3 x 10^ ti pairs, respectively, as given in Table |3.1| . Therefore, the 
LHC would be able to probe the CP asymmetry of the top quark at the level of a few 
percent. A similar number of the ti pairs is required in electron collisions to probe 

the CP asymmetry at the same level. Thus, for a = 500 GeV e~e~^ collider, an 
integrated luminosity of about 10^ — 10^ fb~^ has to be delivered. This luminosity is 
at least a factor of 100 higher than the planned next linear colliders. We note that 
although the initial state in a pp collision (such as at the LHC) is not an eigenstate 
of a CP transformation, these CP-odd observables can still be defined as long as the 
production mechanism is dominated by gg fusion. This is indeed the case for ti pair 
production at the LHC. 

In the SM, the top quark produced via the P^-gluon fusion process is about one 
hundred percent left-hand (longitudinally) polarized, see Appendix A. Given a polar- 
ized top quark, one can use the triple product correlation, as defined in Equation ( |7.2| ), 
to detect CP violation of the top quark. For a polarized top quark, one can either 
use at X ph or p]^^^ x to define the decay plane of t ^ bW{^ £"*"//). Obviously, the 
latter one is easier to implement experimentally. Define the asymmetry to be 

^ _ N{i^ out of the decay plane) — N{£'^ into the decay plane) 
*° N{i+ out of the decay plane) + N{i+ into the decay plane) 

If Aio is not zero, then the time-reversal T is not conserved, therefore CP is violated 
for a CPT invariant theory. Due to the missing momentum of the neutrino from the 
decay of the VT-boson, it is difficult to reconstruct the azimuthal angle {4>w) of the 
W-boson from the decay of the top quark. Once the angle (pw is integrated over, 
the transverse polarization of the top quark averages out and only the longitudinal 
polarization of the top quark contributes to the asymmetry Aio- Thus, the asymmetry 
Aio can be used to study the effects of CP violation in the top quark, which in the 
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SM is about one hundred percent left-hand (longitudinally) polarized as produced 
from the W^-gluon fusion process. To apply the CP-violating observable Aio, one 
needs to reconstruct the directions of both the t and b quarks. It has been shown in 
Reference that it takes about 10^ — 10^ single-top events to detect CP violation 



at the order of ~ 10 ^ — 10 ^. 

For rrit = 180 GeV at the Tevatron, the Di-TeV and the LHC, an integrated 
luminosity of 10, 100 and 100 fb"^ will produce about 2 x 10^ 1.4 x 10^ and 2 x 10^ 
single-t or single-t events, respectively, from Table |3.1|. At the NLC, the single top 



quark production rate is much smaller. For a 2 TeV electron collider, the cross sections 
for e~e~^ e'vjh and e"'"7 Ugtb are 8 fb and 60 fb, respectively [^. Hence, it 



will be extremely difficult to detect CP violation effects at the order of < 10 ^ in the 
single-top events produced in electron collisions. 

A few comments are in order. First, to extract the genuine CP- violating effects, 
we need to study the difference in the asymmetry Aio measured in the single-t and 
single-t events because the time-reversal violation in Aio of the t (or t) alone could be 
generated by final state interactions without CP-violating interactions. Second, the 
detection efficiency for this method is not close to one, so a good understanding of 
the kinematics of the decay products and how the detector works are needed to make 
this method useful. 

The asymmetry Aio belongs to the second category of CP-violating observables 

and is CP-odd and CPT-even. Consider another asymmetry At which belongs to the 

first category of CP-violating observables and is CP-odd and CPT-odd. Using At for 
detecting CP-violating effects is to make use of the fact that pp is a CP eigenstate; 
therefore, the difference in the production rates for pp — >■ tX and pp — >■ iX is a signal 
of CP violation. This asymmetry is defined to be 

. ^ a(pp ^ tX) - ff(pp ^ iX) 
~ cr{pp ^ tX) + a{pp tX) ■ 

As discussed in Chapter 6, the production rate of pp — > tX is proportional to the 
decay rate of t — *• bW~^ and the rate of pp iX is proportional to the rate of 
t — > bW~ . This imphes that At = Atw, cf. Equation (|73| ). There have been quite 
a few models studied in the literature about the asymmetry in Abw- For instance, 
in the Supersymmetric Standard Model where a CP-violating phase may occur in 
the left-handed and right-handed top-squark, Abw can be as large as a few percent 



depending on the details of the parameters in the model [^] . 

Next, let's examine how many top quark events are needed to detect a few percent 
effect in the CP- violating asymmetry At- Consider t — >■ bW~^ bi'^u, where i = 
COT fi. Define the branching ratio as the product of Br(t —>■ bW~^) and Br(iy+ 
£+z/), where Br{W^ £+z/) is 2/9. (Br(t — > bW^) depends on the details of a 
model and is almost 1 in the SM.) Let us assume that the efficiency of 6-tagging 
(^btag) is about 30% and the kinematic acceptance (e^) of reconstructing the single- 
top event, pp tX — > bW~^X bi'^uX, is about 50%. (See, a Monte Carlo study 
in Chapter 8.) The number of single-t and single-t events needed to measure At is 
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Thus, to measure At of a few percent, A/^ has to be as large as ~ 10^, which corre- 
sponds to an integrated luminosity of 100 fb~^ at the Di-TeV. 



Chapter 8 

A Monte Carlo Study 



It was shown in Reference p4| that due to the characteristic features of the transverse 
momentum and rapidity distributions of the spectator quark which emitted the virtual 
W an almost perfect efficiency for "kinematic b tagging" can be achieved. In addition, 
the ability of performing 6-tagging using a vertex detector increases the detection 
efficiency of a heavy top quark produced via the W^-gluon fusion process. In this 
chapter we show that this process is useful at the Tevatron with the Main Injector. 

We also estimate results for = 4 TeV at the Di-TeV and for = 14 TeV at the 
LHC in separate subsections. 

To show that a heavy top quark produced from the TV-gluon fusion process can 
be detected at the Tevatron, we performed a Monte Carlo study on the W + 2 jets 
mode of the (2 — 2) process 

qb qt{^ bW^{^ tv)) (8.1) 

with = e+or/i+. More specifically, we assume that the 6-quark jet from the top 
quark decay can be tagged so that the decay mode of interest is identified to be 
W + b + jet. Throughout this study, we assume that the efficiency of the 6-quark 
tagging is 30% for > 30 GeV with no misidentifications of a 6-jet from other 
QCD jets. For clarity we only give rates for top quark (not including top-antiquark) 
production in this chapter, unless specified otherwise. To include the top-antiquark 
production one can refer to Chapter 3 for its production rate as compared with that 
of a top quark. 

For simplicity we only consider the intrinsic backgrounds {i.e., those present at 
the parton level) for the W + b + jet final state, and will not invoke any detailed 
study on effects due to hadronization of the partons or the imperfectness of the 
detectors used in experiments. The intrinsic backgrounds in the SM for the mode 
W + b + jet are the electroweak-QCD process q'q ^ W + b + b and QCD process 
qq, gg ti —>■ W + b + jet. We will show that the dominant backgrounds for the 
single-top signal come from the electroweak-QCD processes (as shown in Figure |8 . 1| ) 

ud, OS ^bbW+{^ tu). (8.2) 

The other backgrounds such as eg —>■ bW^ are suppressed due to the small CKM 
matrix element \Vcb\ — 0.03 to 0.048 0. As done in the previous sections, we will give 
our numerical results in this section for either a 140 GeV or a 180 GeV top quark. 
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8.1 Tevatron with VS' = 2 TeV 



At the Tevatron (pp, = 2 TeV) the single-top production rate from the VT-gluon 
fusion process is 1(2) pb, as shown in Table For the final state of Equation ( |8.1| ), 
the branching ratio of t — > bW^{^ i'^u) for = e"*" or /i+ is about 2/9 in the SM. As 
discussed in Chapter 3, we found that after properly treating the b quark as a parton 
inside the proton (or antiproton) the total rate for the H^-gluon fusion process is 
about 30% smaller than that of the (2 2) process. Hereafter, we shall rescale all 
the numerical results of our analysis for the (2 2) process to the total event rate 
of the VT-gluon fusion process by multiplying them by a factor of 0.7. 

To show that a 180 (140) GeV top quark produced from this process can be 
detected at the Tevatron, we first impose the following kinematic cuts: 



P| > 15 GeV, 

> 15 GeV, 

> 35 GeV, 



|r/«| < 3.5, 
\V^\ < 2, 
\v'\ < 2, 



> 15 GeV, ARgb > 0.7 . 



(8.3) 



The efficiency of these cuts for a 180 (140) GeV top quark is 32%(53%). Including the 
reduction factor from the assumed 6-tagging efficiency, 30%, the signal rate is found to 
be about 0.045 (0.063) pb. In Equation ( ^.3|) , Pt stands for transverse momentum, rj 



for pseudo-rapidity, fix for missing transverse momentum, and AR = ^ {SrjY + (<50)^ 



with being the azimuthal angle. It is important to note in Figure |8.2| that the typical 
rapidity of the spectator jet in the signal event is about 1.6 although almost all the 
signal events have \ri'^\ < 3.5 The distribution of rj'^ is asymmetric because the 
Tevatron is a pp collider. To produce a heavy top quark, which decays to a positively 
charged lepton, the valence quark from the proton is most important, implying a large 
probability for t]'^ to be positive. (We define the positive z-direction to be the proton 
moving direction in the laboratory frame.) Similarly, a top-antiquark produced from 
the ly-gluon fusion process would prefer a negative rj'^ due to the large up-antiquark 
PDF inside the antiproton. 

In the W + bb background process, the bb pair comes from a virtual gluon con- 
version, therefore its rate is highly suppressed if the invariant mass of the bb pair is 



u, c 



d. s 




u, c 



d. s 




Figure 8.1: Diagrams for ud, cs —>■ bbW 



46 



15.0 



12.5 



10.0 



7.5 



5.0 



2.5 



1 1 1 r 



0.0 




-4 



Figure 8.2: The rapidity distribution of the spectator quark g, after cuts in Equa- 
tion ( ^.31) , for the signal q'b — > > bW^{-^ ^'^^ of the spectator quark h 
for the major background q'q bbW~^{^ i'^u) (dots), for = 180 GeV (sohd) and 
140 GeV (dash), at the Tevatron. (The vertical scale is arbitrary, but the relative size 
among these curves are absolute.) 
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large. Since both b and b have about the same transverse momentum (Pt) in the 
background event, the requirement of > 35 GeV effectively forces a similar Pt 
cut on b. This generates a large invariant mass of b and b {i.e., the virtuality of the 
gluon), strongly suppressing the background rate. In contrast, in the signal event the 
final parton q (from q', after emitting a virtual W) typically has a smaller Pt than 
the 6-quark (from the decay of a heavy top quark). Typically, in the signal event, 
P|. ~ mj/3. In Figure O we show the Pt distribution of the tagged b from t in both 



the signal and the background events. Hence, demanding an asymmetric cut on Pt 
{i.e., Pt > 15 GeV and P^ > 35 GeV) will suppress background effectively and keep 
most of the signal events. This is why an asymmetric cut on Pt was used in our 
analysis to suppress the major background process W + bb. To compare the efficiency 
of this asymmetric cut in Pt, we note that using P|- > 15 GeV along with all the 
other cuts in (|8.3|) yields a signal-to-background ratio (S/B) of about 1/3(2/3). Re- 
quiring > 35 GeV excludes about 60% of the background events sacrificing about 
10%(30%) of the signal. 

After imposing the kinematic cuts in Equation ( ^.3|) , we found that S/B ^ 0.9(1.3). 
However, the signal-to-background ratio can be further improved by imposing 

cos^^, >-0.4. (8.4) 

Because the top quark produced from the VT-gluon fusion process is left-hand polar- 
ized, tends to move against the moving direction of the top quark in the center- 



of-mass frame of q and t, cf. Figure |6.1| . However, in the background event, the 



distribution of cos^^^, as shown in Figure |8.4| , is almost fiat after imposing the cuts 



of ( p.3|) . {dig = n — Oi, where 6i is the the polar angle of in the rest frame of t 
defined in the center-of-mass frame of q and t.) 

To calculate cos^^g, the P^ information must be constructed. Since both and 
u come from a real boson, we can use the VT-boson mass constraint 



M^ = {pe + p,y (8.5) 

and the ^t information to specify the longitudinal momentum P^ of the neutrino. 
There are two solutions for P^ and typically, both of them are physical solutions for 
a signal event. Therefore, one has to fix a prescription to choose the one which will 
most likely give the correct distribution of the invariant mass of £+,z/ and b. We 
choose the solution which has the smaller |P^|. Here we exploit the fact that the W 
has finite width. If a physical solution for P^ is not found with = 80 GeV, we 
generate a resonant mass of the ly-boson using a Breit-Wigner distribution. We use 
a full half-maximum width of the ly-boson, where = 2.1 GeV, and solve for P^, 
repeating the algorithm for up to three trials if necessary. We found that the survival 
probability for finding a solution using this algorithm is about 90%, and the difference 
between this solution and the value actually generated by the Monte Carlo generator 
is a Gaussian distribution peaking at with a width about the order of Tw After 
the additional cut imposed on cos9eq, we obtained S/B ~ 1.2(1.8). About 55%(40%) 
of the total signal event rate remains after applying the cuts ( p.3| ) and ( p.4|) to the 

process ( |8.1| ). We conclude that for an integrated luminosity of Ifb^^ at a 2TeV pp 



collider, there will be about 75(105) signal events detected with a significance S/vB 
of about 10(14), including both the single-t and single-t events as defined in ( ^.1|) . 
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Figure 8.3: Pt distribution of the h quark, after requiring P|. > 15GeV along with 
all the other cuts in (p.3|), for the signal q'h qt{-^ hW^{-^ and the major 

background q'q hhW^{—^ ^'^^), at the Tevatron. 
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Figure 8.4: The cos 9iq distribution prior to cut Equation (|8.4|) for the signal q'b 
qt{-^ hW^{-^ i'^i')) and the major background q'q — >■ bbW~^{^ ^^^), at the Tevatron. 
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Figure 8.5: The rrit distribution after the cuts Equation ( |8.3| ) and Equation ( |8^ ) for 
rrit = 180 GeV (sohd) and 140 GeV (dash) at the Tevatron including both the signal 
and background events with — > or/x^*". 
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To measure the mass of the top quark, we calculate the reconstructed invariant 
mass {mt) of the top quark using 



= {pb + Pe + Puf- (8.6) 

The distribution of rrit, including both the signal and the backgrounds, is shown in 
Figure |8.5| , in which a clear mass peak appears unmistakably. Therefore, we conclude 
that the top quark can be detected and studied via this process at the Tevatron. In 
Figure |575|, other less important backgrounds, such as ti events, were also included. 



The tt background is not important after vetoing the events with more than 2 jets 
2^ . To support this we did a study for the two decay modes of tt: 

ti hW+ (^ tv) bW^ (^ q'q) (8.7) 

and 

ti^bW^i-^ i^u)bW-{^ i^u), £2 = eor/i. (8.8) 
For both modes we require one of the jets to be a 6 or b. 

Consider first the i —>■ bq'q decay mode. We require at least one of the q'- and 
g-jets and b (if b is tagged) or b (if b is tagged) to be within rapidity 3.5, otherwise 
we reject the event. If any one (and only one) of the untagged jets is within rapidity 
3.5 we call it the spectator jet and then apply our cuts. If two of the untagged jets 
are within rapidity 3.5, then we require their AR separation to be less than 0.7 to 
classify them as one spectator jet. When all three untagged jets are within rapidity 
3.5 we choose the jet with the largest Pt and check its Ai? separation with the other 
two. If the lower Pt jets are within Ai? = 0.7 of the high Pt jet, we call this the 
spectator jet and apply our cuts, otherwise, reject the event. Recall from Tabl e ^7\\ 
for a 180 GeV top quark, the ti rate is about 4.5 pb. After applying the cuts (|8.3| ) 
and (^.4|) and including the branching ratio for this mode, || = ^5 the event rate, of 

approximately 3 x 10~^ pb, is very small as compared with the signal rate of 0.075 
pb (including t and i). This is because for most of the i — >• bq'q decay modes all three 
jets are within rapidity 3.5 and Ai^^j and Ai^g^ are in general large . 

For the i decay mode we require that £2 be undetected. Specifically, if 

£2 is within rapidity 2 with P^ > 15 GeV we reject the event. If 2 < \ri^^\ < 3.5 we 

require P^^ to be less than the minimum Pt for detecting leptons, i.e., < 15 GeV in 
accordance with Equation ( ^.3| ). After the cuts in ( ^.3[ ), this mode (with branching 
ratio If = ^) already suffers, being about 5% of the signal rate. This mode suffers 
another factor of about two loss to the failure of reconstructing P^ due to the presence 
of two neutrinos in the final state. After imposing the cos Oiq cut the rate for ti —>■ 

bW~^{^ £+z/) bW~{^ ^2'^) is about 3 x 10~^ pb which is about a factor of 25 smaller 
than the signal rate. Hence the dominant background (of the same order as the signal 
rate) comes from the electroweak-QCD processes as given in Equation ( |8.2| ). 

As summarized in Figure even with the very minimum kinematic cuts of (|8.^ 



and ( ^.4|) the single-top signal can already be detected, assuming a perfect detector 



with 6-tagging efficiency of 30%. To incorporate the effects of detector efficiencies. 
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we smear the final state parton momenta using a Gaussian distribution with 



{AE/E)i = 15%/ Ve, and {AE/E)g^b = 50%/Ve. (8.9) 



The nit distribution becomes shghtly broader as shown in Figure |8.6| ; however, both 
the signal and the background rates are almost the same as those obtained with 
a perfect detector. As discussed in the previous chapters, the distribution of the 
invariant mass rriM is extremely useful in either determining the mass of the top quark 
or measuring the form factors of t-b-W . For completeness we show in Figure |8.7| the 
distribution of mM with or without smearing, after all the above analysis. Because 
the 6-jet is required to have large Pt, cf. Equation (^.3|) , so AE/E for the 6-jet 
is small, therefore the two mbi distributions do not differ much and the difference 
becomes smaller for larger m^^. Thus, the position of the bump at large rrihi in 
the signal events remain a good signature for detecting the single-top signals and 
determining the mass of the top quark (discussed in Chapter 4) or the couplings of 
t-b-W (discussed in Chapter 6). 

We note that the data sample obtained after all the above analysis can be further 
purified at the cost of somewhat reducing the signal rates. This can be easily done, 
for instance, by noting the distinct differences between the signal events and the 
background events from the distributions of rapidity of the spectator jet (Figure |8^) , 
transverse momentum of the bottom quark (Figure p.3|) , and the angular correlation 
cos 9iq (Figure p^ ) due to the polarization of the top quark in signal events. However, 
for a more realistic simulation, one should also consider the possibility of having a 
charm-jet (or even an ordinary QCD-jet) faking a bottom jet in 6-tagging so that the 
actual background rate measured by the detector would be larger than that given 
here. This is outside the scope of our parton level study. 

As discussed in Chapter 3, another process which produces a single-top is the W* 
production 

q'q ^W* ^ tb. (8.10) 
At 2 TeV for a 180 (140) GeV top quark, the W* production rate is about 1/5(1/3) 



of the ly-gluon fusion rate. Applying the kinematic cuts defined in Equation 
we find that the W* process passes with about the same efficiency as the ly-gluon 
fusion process. However, there are a few obvious differences in the kinematics of their 
final state partons. First, in the W* event, there are two 6-jets (one for b from t 
decay and another for b from production), therefore there is a 50% chance of tagging 
the wrong b and giving the wrong reconstructed top quark invariant mass, as defined 
in Equation (|8.6| ). To improve the invariant mass distribution of the top quark, one 
has to be able to distinguish a 6-jet from a 6-jet by making further selections at the 
cost of reducing the single-top rate from this process. (Some of the techniques have 
been discussed in Chapter 4.) Second, the rapidity distribution of the spectator jet 
{i.e., 6-jet if 6-jet identified) in the W* event peaks around zero {i.e., central, as 
shown in Figure |8.8| ) because tb is produced through the s-channel process.^ This is 
in contrast to that in the py-gluon fusion event where the rapidity distribution {rj'^) 

^ We note that the rapidity distribution of the 6-quark in the W* tb event is shghtly asymmetric 
around zero. It favors a shghtly negative rapidity. (Recall that in the M^-gluon events for producing 
single-f, the rapidity of the spectator quark q favors positive values.) This is similar to the lepton 
rapidity asymmetry observed in the pp t^v events due to the ratio of the down-quark 

and the up-quark parton distributions inside the proton and the anti-proton. 
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Figure 8.6: Same as Figure |8]^, but with detector resolution effects as described in 
Equation ^]9[ 
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Figure 8.7: rriM distributions without (top) or with (bottom) smearing, after all the 
above analysis for the Tevatron. 
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Figure 8.8: The rapidity distribution, after Equation (|8.3| ), of the spectator jet {i.e., 
6-jet if 6-jet identified) in the W* event for mj=180 GeV (sohd) and 140 GeV (dash) 
against Figure |8]^ (dots), at the Tevatron. 
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of the spectator jet (labeled as g-jet) is asymmetric and less likely to be around zero. 
Third, the polarization of the top quark produced from the W* process is not purely 
left-hand polarized as in the case of the q'b qt process. For a 180 (140) GeV top 
quark, the ratio of the event rates for producing a left-handed top versus a right- 
handed top in the W* event is about 3.5(3.4). In Figure |8l9| , we show the production 



rate for a left-handed or a right-handed top quark from the W* process and compare 
them with that from the H^-gluon fusion process as a function of m^. Because the 
top quark is not 100% polarized in the W* process, the angular correlation of 
and the spectator jet will not be as strong as that in the W^-gluon fusion process. In 
Figure |8.10| we show the distribution of cos O^q in W* events for a 140 and 180 GeV top 



quark. (Here, g-jet denotes the spectator jet.) Following through the previous analysis 
done for the W^-gluon fusion events we found that W* production compliments the 
VT-gluon fusion process by increasing the single-top production rate by about 1O%.0 
Therefore, its contribution to our final results of various distributions is small. 

In conclusion we found that at the Tevatron (pp, \fS = 2 TeV) the single-top 
production rate from the H^-gluon fusion process after including the branching ratio 
for t hW+{-^ i+u) is about 0.22(0.44) pb for a 180 (140) GeV top quark, where 

= e'^OTfi^. The kinematic acceptance after the kinematic cuts (|8.3|) and (|8.4| ) is 
about 55%(40%). Assuming a 30% 6-tagging efficiency we concluded that the single- 
top event rate from the VT-gluon fusion process is about 0.036(0.052) pb. For an 

integrated luminosity of 1 fb~^, this yields 36(52) reconstructed single-top events. 
(To include top-antiquark production, a factor 2 should be included.) The dominant 
background process is the electroweak-QCD process W + bb whose rate is about 
80%(60%) of the signal rate in the end of the analysis. The ti events are not as 
important to our study. The ti rate for a 180 GeV top quark is only 0.4% and 4% of 
the signal rate for its lepton-|-jet and di-lepton mode, respectively. In both Figs. ^]5| 
and ^]6|for the distribution of rrit, we have also included another single-top production 
process (a single-top produced from W*) which increases the single-top rate by about 
10%. 

Let us make a side remark about the dominant background q'q W + bb before 
we close this section. In the above analysis we did not include the possibility of having 
an additional QCD jet from either the radiation or the conversion of the incoming 
quark jet ( q' or q). The concern is that this jet may be identified as a forward jet 
which would fake the single-top signal event. In this case, b and b in the W + bb+ jet 
background event have to both fall into a cone of AR = 0.7 in order to fake the tagged 
6-jet (only one b) in the single-t event (cf. Equation (p.l|) ). To examine the possibility 
for this to happen, we have applied the eikonal approximation |^0[ to calculate the 
rate ofW + bb+ jet from the square of the W + bb amplitude. (The results are shown 
in Appendix E.) 

After the basic kinematic cuts: 

P^"* > 15GeV, Ir]^"'] < 3.5, 

> 15GeV, \r]^\ < 2, 

Pf^ > 15GeV, \7]'''^\ < 3.5, (8.11) 



^ The W* production rate is about one fifth of the M^-gluon fusion rate for a 180 GeV top quark, 
and the kinematic acceptance of the W* event is about half of that of the VF-gluon fusion event. 
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the rate for W + bb+ jet is already about a factor of 5 smaller than that for W + bb. 
Naively, one might expect a factor of as{Mw) (~ 0.1) suppression factor for emitting 
an additional QCD jet (quark or gluon) in the hard scattering process. However, as 
compared to the large invariant mass (M^y^^) of the W + bb system a jet Pt of 15 GeV 
may be small enough to generate large logs, such as \ia{M-iYi,i/ Pt), in the amphtudes. 
Hence, because ln(M^^5/PT) is not negligible, the rate for W + bb + jet is not 
suppressed by a factor of 10 relative to the rate of + bb, but only a factor of 5. 



To see how often b and b will fall into a cone of AR = 0.7 we show in Figs. |8.11 



and |8.12| the Ai?^^ and the M^^ d istributions in W + bb + jet events^ after applying 
the kinematic cuts listed in ( ^.11| ). The same distribution in W + bb events is also 
shown for comparison. The ARf^i distributions look alike, and the Mj,^ distribution 
falls slowly as Mg, increases. Also, the W + bb+ jet event prefers a larger M^j because 
the Pt of b and b are larger in this process than that in the W + bb process. We find 
that only about 20% of the W + bb + jet events can possibly fake the single-t event 
by having b and b inside the same jet-cone and thus fake a tagged 6-jet. Hence, the 
additional background rate from W + bb+ jet events is about a factor of | x 20% = 4% 
of the electroweak-QCD background rate. Although our estimate is not precise, we 
believe our conclusion for this additional background should hold within a factor of 
2. Hence, this additional background is negligible at the Tevatron. However, it can 
be important at the LHC Because the energy of the LHC collider is much higher, it 
is more likely to have additional radiation in the event and to boost the bb system to 
make them closer and thus fall into the same jet-cone. 



8.2 Tevatron with = 4 TeV 

Here we present our results for a possible upgrade of the Tevatron with = 4 TeV. 
After the following kinematic cuts: 

P|>15GeV, |r/«|<3.5, 

P|>15GeV, \r]^\<2, 

P^>30GeV, |r/''|<2, 

^T > 15 GeV, APgfc > 0.7 , (8.12) 
the signal rate is about 0.28(0.37) pb. (The efficiency of these cuts is 45%(56%).) 



In Figure |8.13| the typical rapidity of the spectator jet in the signal event is about 
2, and almost all the signal events have \ri'^\ < 3.5. An asymmetric cut on Pt was 
used once again to suppress the major background process W + b + jet. Demanding 
Pt > 15 GeV along with the other cuts in Equation (p.l2|) , the signal-to-background 
ratio (S/B) is about 1.1(1.7). We show in Figure ^.14| the Pt distribution of the 
tagged b from t. Requiring P^ > 30 GeV excludes about half of the background events 
sacrificing about 6%(20%) of the signal. After all the cuts listed in Equation (|8.12| ), 
S/B ~ 2.3(3.0) with the signal rate at 0.28(0.37) pb. 

We show in Figure |8.15| the distribution of cosOgq. After applying the cut 



the ratio S/B ~ 2.9(3.8) with the signal rate of 0.22(0.29) pb. In the end of the 
analysis there are about 2200(2900) single-t events for an integrated luminosity of 10 
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fb~^ at = 4 TeV (a pp collider) with a significance S/v^ of about 80(105). The 
kinematic acceptance of the signal event is about 43%(34%). Note that in all the 
above rates we have included the reduction factor from a 30% 6-tagging efficiency. 
In Figure ^.16| , we show the reconstructed invariant mass (m^) of the top quark. 
Once again, to incorporate the effects of detector efficiencies, we smear the final state 
parton momenta as in ( p.9| ). The distribution becomes slightly broader as shown 
in Figure however, both the signal and the background rates are almost the 

same as those obtained with a perfect detector. 



8.3 LHC with = 14 TeV 



Here we present our results for the LHC with ^/S = 14 TeV. After the following 
kinematic cuts 

P|>40GeV, 1<|77«|<4, 

P-^>40GeV, \v^\<2, 

P^>40GeV, |r/''|<2, 

> 40 GeV, ARgb > 0.7 , (8.13) 

the signal rate is about 0.44(0.24) pb. (The efficiency of the cuts is 94%(98%). We 
still assume a 30% efficiency for the 6-tagging at the LHC.) Here we did not impose 
a smaller cut because a lower Pt jet will be more difficult to be identified at 
the LHC. (A typical QCD event at the supercollider will be engulfed by soft gluon 
radiation.) Since the signal event yield is large at the LHC, we decided to purify 
our data simply by requiring a large P^ cut. Notice that the rate for rrit = 180 GeV 
is larger at the LHC than that for = 140 GeV after our cuts, opposite to the 
behavior of the rate at 2 TeV and 4 TeV. This is due in part to less sensitivity to 
at higher energies, but mainly because the b from top decay is much harder for larger 
top mass and thus is less sensitive to the cut of P^ > 40 GeV. The typical P^ for 
= 180 (140) GeV is 60(40) GeV. We show in Figure (PTT^ ) the Pt distribution of the 
tagged b from t. 



In Figure |8.19| the typical rapidity of the spectator jet in the signal event is about 
3, but a cut on \ri'^\ < 4 keeps almost all the signal events. Excluding the |?7^| > 1 
cut in (|8.13|) the signal-to-background ratio (S/B) is about 10(7). Requiring \ri'^\ > 1 
excludes about 40% of the background events sacrificing about 15%(20%) of the 



signal. After the kinematic cuts in (^.13|) , the ratio S/B ~ 25(14) 



We show in Figure g.20 the distribution of cos Ogg at the LHC. After the cos 9ig cut 



S/B ~ 40(20). In the end of the analysis, there are about 30,000(15,000) single-top 
events for an integrated luminosity of 100 fb~^ at ^/S = 14 TeV (a pp collider) with 



a significance S/vB of about 32(16). Hence, about 4%(2%) of the total signal event 
rate remains. This is thus the kinematic acceptance for the signal process (|8.1|) . In 
Figure ^.21] , we show the reconstructed invariant mass {rrit) of the top quark for a 



perfect detector. Once again, to incorporate the effects of detector efficiencies, we 
smear the final state parton momenta as in (|8.9|) and reconstruct the rrit distribution 
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for the LHC is shown in Figure 8.22| . 

Notice that the above analysis for the LHC is less reliable because the energy of 
the collider is much higher and therefore it is more likely to have additional soft- 
jets accompanying the signal and the background events. As discussed at the end of 
section 7.1, it would be more reliable to use a full event generator such as ISAJET 
pl| , PYTHIA |]62[ or HERWIG |Q for this study because these generators contain 



radiation from either the initial or final states. However, these generators currently do 
not have the correct angular correlations in £ and jets, as discussed in this analysis. It 
would therefore be important in the future to improve these generators to incorporate 
the polarization effects of the top quark and the ly-boson for studying physics of the 
top quark in hadron collisions. 



60 




I \ \ \ \ I \ \ \ \ I \ \ \ \ I I I I I I 

100 150 200 250 300 

Top Mass mt(GeV) 



Figure 8.9: The production rate for a left-handed (long dash) or a right-handed 
(short dash) top quark from the W* process. The upper solid line is the total rate 
for the H^-gluon fusion process, the lower solid line for the W* process. 
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Figure 8.10: The distribution of cos^^g in W* event for mt=180 GeV (solid) and 140 
GeV (dash) against Figure (dots), at the Tevatron. 
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Figure 8.11: Ai?^^ distributions in W + bb + jet (solid) and W + bb (dash) events 
after applying the cuts listed in (|8.11|) . 
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Figure 8.12: distributions in W + bb+ jet (solid) and W + bb (dash) events after 
applying the cuts listed in (|8.11|) . 
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Figure 8.13: The rapidity distribution of the spectator quark g, after cuts in Equa- 
tion (|8.12| ), for the signal q'b qt{-^ hW^{-^ ^ind of the spectator quark h 
for the major background q'q — >• hhW^{-^ i^u) (dots), for = 180 GeV (sohd) and 
140 GeV (dash), at the Di-TeV. 
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Figure 8.14: Pt distribution of the h quark, after requiring P|. > 15 GeV along with 
all the other cuts in ( ^.12| ), for the signal q'h — > qt{—^ hW^{-^ ^^^))) cind the major 
background q'q hhW^{—^ ^^v)-, at the Di-TeV. 
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Figure 8.16: The rrit distribution after the cuts Equation ( |8.12| ) and Equation 
for mt = 180 GeV (sohd) and 140 GeV (dash) at the Di-TeV, including both 
signal and background events with 



or . 
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Figure 8.17: Same as Figure |8.16| but with detector resolution effects as described 
in Equation ( p.9|) . 
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Figure 8.18: Pt distribution of the b quark, after cuts in ( p.l3| ), for the sig 
q'b qt{—* bW^{-^ ^^^))) the major background q'q —* bbW^{—^ at 
LHC. 
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Figure 8.19: The rapidity distribution of the spectator quark g, after cuts in Equa- 
tion (|8.13| ), for the signal q'b qt{-^ hW^{-^ ^ind of the spectator quark h 
for the major background q'q bbW^{^ i^u) (dots), for = 180 GeV (sohd) and 
140 GeV (dash), at the LHC. 
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Figure 8.20: cosO^q distribution for the signal q'h — >• >■ hW^{-^ ^^^)) ^^'^ 
background g'g — >■ hhW^{-^ t^v) at the LHC. 
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Figure 8.21: The rrit distribution after the cuts Equation ( |8.13| ) and Equation ( |8.4| ) 
for rrit = 180 GeV (sohd) and 140 GeV (dash) at the LHC including both the signal 
and background events with 



or fi"^ . 
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Figure 8.22: Same as Figure |8.21| , but with detector resolution effects as described 
in Equation ( p.9|) . 



Chapter 9 

Discussions and Conclusions 



We discussed the physics of top quark production and decay at hadron colhders, such 
as the Tevatron, the Di-TeV and the LHC We showed how to use the invariant mass 
distribution of nihi to measure the mass and the width of the top quark, produced 
from either a single-top or a ti pair process. It has been shown in Reference |^ that 



the distribution of mf,^ is not sensitive to radiative corrections from QCD interactions. 
Thus it can be rehably used to test the polarization of the W^-boson from t decay 
(hence, test the polarization of the top quark from the production mechanism) and 
to measure the mass of the top quark using the observed value of /Long (the fraction 
of longitudinal W^s from top decays). We also discussed how well the couplings of 
t-b-W vertex can be measured to probe new physics, and how well the CP properties 
of the top quark can be tested in electron or hadron colliders. 

In Reference we showed that an almost perfect efficiency for "kinematic b 



tagging" can be achieved due to the characteristic features of ly-gluon fusion events. 
In addition, the ability of performing 5-tagging using a vertex detector increases the 
detection efficiency of a heavy top quark produced via the ly-gluon fusion process. 

A detailed Monte Carlo study on how to identify the characteristic features of the 
signal events {i.e., the transverse momentum and the rapidity distributions of the 
spectator quark which emitted the virtual W) and therefore suppress the background 

events was performed in Chapter 8.QFor an integrated luminosity of 1 fb~^, there will 
be about 75 (105) single-t or single-t events reconstructed in the lepton+jet mode for 

■rrit = 180 (140) GeV at y/S = 2 TeV. (The branching ratio of — > e, or /i is included, 
and the 6-quark tagging efficiency is assumed to be 30% for > 30 GeV with no 
misidentifications of a 6-jet from other QCD jets.) The dominant background process 
is the electroweak-QCD process W + bb whose rate is about 80%(60%) of the signal 
rate in the end of the analysis. The ti events are not as important to our study. The 

results for = 4 TeV at the Di-TeV and for = 14 TeV at the LHC were also 
discussed. 

Although the W* — * tb rate in the SM is not as large as the VT-gluon fusion rate 
for producing a heavy top quark, it remains a complementary process for probing 



^ The fortran code, ONETOP, used for this study is available by request. In Appendix F we 
briefly describe the processes included in this program. 
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new physics in the single-top quark event. The W* process is particularly useful for 
detecting new physics through some possible high mass resonance in the theory. In 
that case, its rate will be highly enhanced by the resonance effects. We however did 
not study such a possibility in this work because its rate depends on the details of 
the models considered. 



Appendix A 

Helicity Amplitude Method 



In this appendix we would like to display the rules for doing calculations at the ampli- 
tude level using the Helicity Amplitude Method. The method breaks down the algebra 
of four-dimensional Dirac spinors and matrices into equivalent two-dimensional ones. 
This algebra is easy to program and more efficient than computing the Dirac algebra 
as it stands. All diagrams are summed and squared numerically. 

In what follows we introduce the Weyl representation of Dirac spinors and matri- 
ces. We also include several example calculations to illustrate the finer points of the 
method. Throughout this paper we use the Bjorken-Drell metric 

g^"" = dtag {1,-1,-1,-1). (A.l) 

The four-momenta have the form in spherical coordinates: 

= {E, \p \ sin 9 COS0, \p \ sin 9 sincp, \p \ cos9) (A. 2) 

with E"^ — \pf = m?. We define the right-hand (R), left-hand (L) and longitudinal 
(0) polarization vectors for a spin-1 field as Q 

e^^-^ = — ^(0, i sin — COS0 cos^, — i COS0 — sin0 cos5, sin 5) 

= --^(0, z sine/) + COS0 cos6', — z cos(/) + sin0 cos6', — sin6') (A.3) 

efn^ = — (|p I, sin6' cos(/), sin6' sin0, _E cos6'). 

The above equations satisfy the identities e^^^^ = —£(1)' ^(o) ~ ^(o)' Pm^(7i) ~ 
^^^^(h') = -^hh', for h, h' = R,L or 0. 

^ For a massless spin-1 field, only the right-handed and the left-handed polarizations are physical. 
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In four component form we define the following. In the Weyl basis Dirac spinors 
have the form 



where for fermions 



and anti-fermions 



with w± — E ± |p |. 



The XA/2's are eigenvectors of the hehcity operator 

h^p- a,p^pl\p\ (A.7) 
with eigenvalue A where A = +1 is for "spin-up" and A = — 1 is for "spin-down". 

/ cos^/2 \ / -e-^<^ sin \ 

^ \ sin QI2 ) ' ^ [ cos ) ' ^^'^^ 

Later it proves useful to represent XA/2's using bra-ket notation where 

|p+> =Xi/2, \p-) =X-i/2- (A.9) 
Gamma matrices in the Weyl basis have the form 

where aj are the Pauli 2x2 spin matrices 

The chirality projection operators are defined by 

P± = l(l±7^). (A.12) 

Notice that P+ projects out the "right-handed" ("left-handed") component of 
the Weyl spinor effectively reducing the algebra from one involving four component 
spinors and matrices to one involving two component spinors and matrices. 
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Figure A.l: Diagram for the (2 — > 2) process ub ^ dt. 



0\ f^+\ ^ ( 

U- 



p_1}) = 

= J) (J o)=(^10) (A.13) 



P.7^-(.„A.. ^°V-^l-f.° i^)-P.(^':^] (A.14) 



In the Weyl basis ^ has the form 

Po + o 
^Po-a-p J ~ 0;~''^V7- 

where 

7^ = (l,Ta). (A.15) 

Products of these 7^'s have the following useful property when Lorentz indices are 
contracted: 

{l+)ij{lti+)ki = {l'^)ij{li,-)ki = '2[SijSki - Su6kj] (A.16) 

and 

{l+)^j{ll,-)kl = {l-)ij{lt,+)ki = 26u6kj, (A. 17) 

where the Roman indices are not vector indices in the usual sense, but are labels 
identifying bras and kets. For instance, for arbitrary kets \i), \j), \k) and |/) we have 

{ithMki = {th^\j){k\^,-\l) = 2{t\l){k\j) = 25uS,j. (A.18) 



Equations ( |A.16| ) and ([A.17|) are simply the two-dimensional version of the well known 
Fiertz identities. 



A.l Helicity Amplitudes for ub dt 

To illustrate the use of helicity amplitudes we calculate the matrix element for the 
(2 2) process ub — > dt which contributes to the total rate for W^-gluon fusion. 
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Figure |A.1| shows the Feynman diagram for this process with the t decay included. 



In this example we use the generalized t-b-W coupling 

= (1 + kl^^) yp- + iiR^^YP+ (A. 19) 

where kl*"*" and kr"'" parameterize deviations from the Standard Model in which 
kl'-^'-^ = and kr'"'" = 0. We calculate the matrix element in the 't Hooft-Feynman 
gauge ignoring for simplicity the factors due to vertices and propagators. We obtain 

M = [u{t)r^'u{b)][u{d)-f^PMu)] (A.20) 

where u, b, d and t are the momenta of the external legs and we retain only the 
mass of the top quark. We use the algebraic properties of the projection operators 
P| = P±,P±P:f = and P±7^ = ■y^P:^ to project out the chirality states. In this 
example the amplitude contains both the left-handed and right-handed currents. For 
pure vector or axial vector currents one must first insert 1 = P+ + P_ or 75 = P-f — P_ 
respectively. Therefore Ai takes the form 

M = {1 + «:^^^)[nUt)7'tn_(6)][nUrf)7^+«-H] 

+ [ul (t)7^! n+ (6)] [ul {d)y,+u^ (n)] (A.21) 
According to Equations ( |A.5| ) and (|A.9| ) we see that 



= J2Eu\u-) 



u_{d) = J2Ed\d-) 



u^h) = j2E,\h-) 



u+{b) = pE,\h+) 



Et-\t\\t+) 



Et + \t\\t-) 



u^t) = { (A.22) 
'Et-\t\\t-) 



Therefore, 



M{+) = (1 + KL'''')^Et - \t\ (£ + u -) 



+ Kif^^Et + \t\ (t + |7'^|6+)(rf-|7^+|M-) 



A^(-) = (l + Ki^^)^E, + |r|(t-|7'^|6-)(t;-|7^+|M-) 

+ /Cfi^^^Et - {i-\y'l\h+){d- 17^+ \u-) (A.23) 
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where A1(±) denotes the amphtude with t hehcity = ±1 and we ignore for now 



the common factor of v^ZE^y^ZEbV^ZErf. 

One is now tempted to move on ahead and contract the Lorentz indices as in 
Equation ([A.16| ). However, contracting 7^'s with the same chirahty introduces extra 



terms into the matrix element. It would be more useful with processes containing 
many branchings from decay, such as in Supersymmetry, if there were a way to utilize 
Equation (A. 17) instead. To this end, we digress a moment. 



Consider some current of the form 

^/±[7MiT---7MnT]V^*± (A.24) 

where n is odd. We note in passing, recalling the properties of the projection op- 
erators, consecutive matrices of the same sign would give zero current. Since Equa- 
tion ( |A.24| ) is just a number, it is identical to taking its transpose as 

^lbl^---ll^Wf±. (A.25) 

We now utilize the following algebraic properties of the Pauli matrices and in partic- 
ular (T2: 

(72(72 = 1, crj = -(72, Cr27jpa2 = 7^± (A.26) 

and define 

= ia2ijl ■ (A.27) 

By inserting pairs of (72 between each pair of objects in Equation ( A.25 ) we therefore 
obtain 

^/±[7mi=f ■ ■ • lf^uT]^i± = ^i±bf.„± ■ ■ ■ 7mi±]V'/± • (A.28) 

It is easy to show that {io'24'i±)~^ = '4'i±' ■ addition, for an even number n of gamma 
matrices we have 

V^/±[7mit ■ ■ ■7/.„±]V^iT = ^iT^b^nT ■ ■ ■7w±]V^/± • (A-29) 

An important result has occurred, which allows us to take advantage of Equa- 
tion ( |A.17| ) avoiding the number of terms that would otherwise occur. For a 



fermion or anti-fermion with momentum p and helicity A, ip± is proportional to either 
or |p— ). It is easy to show that 



\P+) = ~\P-), 

\P-) = 

{p + \ = -{p-\, 

{P^\ = +{p + \- (A.30) 



ia2 acts as a kind of charge conjugation operator on the chirahty states of the Weyl spinors. 
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Finally, recalling Equation ( |A.23| ) the term 

(rf-|7^+|M-) = (M+|7^,_|d+) (A.31) 
riving, via Equation (|A.17|) , 



M{+) = 2(1 + KL^^)\jEt - \t\ {i+\d+){u + \h-) 
+ 2KR^^jEt + \t \ {i + \u-){d-\h+) 



M{-) = 2{1 + KL^^)^jEt + \t \ {i-\d+){u+\b-) 

+ 2KR^^^Et - \t \ {i-\u-){d-\h+). (A.32) 

We now have the matrix element in the form we require for our Monte Carlo package 
ONETOP^, remembering to include coupling constants, propagators, color factors 
and v^v^v^- 

A. 2 Helicity Amplitudes for uh dt in the CMS. 

To illustrate our claim that in the SM (i.e. t^L"^ = and t^n"'" = 0) only the 
left-handed top quark is produced from the ub — > dt process in the d-t center of mass 
frame (CMS), we evaluate the matrix element in terms of CMS variables. Define the 
four-momenta: 

= (71/2,0,0,-71/2) 



= (Vs/2,0,0, Vs/2) 
# = (t, -t sine, 0,-t cos ^) 

= {Et,t sine, 0,t cose) (A.33) 

where Et = {s + m^)/2-\/I, t = {s — mf)/2^/I and we have chosen = to be the 
scattering plane. 

Using these, we obtain from Equations ( |A.5|) , ( |A.8|) and (|A.9| ) 



u+) 

U-) 





-1 

1 




A FORTRAN code. 
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\b + 
\b- 
d + 
d- 
\i + 
\i- 



Therefore, 



{u+\h-) 
{d-\b+) 
{i+\d+) 

{t+\u-) 
{i-\d+) 

{i-\u-) 



Including the common factor 



1 



1 

sin e/2 

-cos e/2 

cos 0/2 \ 
sin e/2) 

cos 0/2 \ 
sin e/2) 

-sin 0/2 
cos 0/2 



-1 

cos 0/2 


cos 0/2 
-1 

-sin 0/2. 



2Ej2EJ2Ea = V - ml) 



and 



we see that 
M{-) 



\lEt + \t\ = 
\lEt - \t\ = 

2 {kr^^) ^Js (s-mf) cos^ 0/2, 



mr 



{AM) 



(A.35) 



(A.36) 



(A.37) 



2(1 + kl^^)Js{s - ml) - 2{KR^^)Jm1{s - mf) sin 0/2 cos 0/2 . 
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u d u d 




Figure A. 2: Diagrams for ug — > dt{-^ hW^{-^ u^)) b . 



Notice in the SM, the top quark is 100% left-hand polarized in the CMS. 

Having outlined the general procedure for calculating amplitudes using the helicity 
amplitude method, we list the matrix elements contributing to single top production 
and the major background Wbb in the Standard Model. We include the decay of 
t — > bW~^ and — > i'^i'e in the final form. 



A. 3 Helicity Amplitudes for ug ^ dt{-^ bW^{- 
t uA) b 



In this and the following sections, we give the diagrams for the process listed, indi- 
cating the rnomentum flow and particle momentum labels: w's are for bosons, 
b's for 6 or 6 quarks, e for e"*", n for z/g and u, d, t and g are for u,d,t quarks and 
gluon, respectively. 

In Figure |A^ , we show the diagrams ioi ug dt{—^ b W^{-^ z/^) ) b. Aii{hg, A^J 
will represent the ith diagram {i = 1,2 from left to right in Figure |A.2| ) where hg 
represents the two transverse gluon polarizations and A^^ represents the two helicity 
states of b. The matrix element in the helicity amplitude formalism for this process 
is 

Mi{hg,T) = 4(62-|^+)(e + | ^- \d+)x 



±yEh,T \bi \{u+\ 4^ Ahg)+ I hi±) Tmb^jEb,± \bi \{u+\ A^g)- \ hi±) 

(g2 - ml) 



M2{hg,T) = ±4V^feiT 1^1 I {b2-\n+){u+\bi±)x 
(e + l f- Ahg)+ f-\d+) +ml{e+\ Ahg).\d+) 



fj,2 _ ^2 



(A.38) 



where we have indicat^ed the four different helicity states involved in this process. We 
keep the mass of the b parton to avoid the case where the b propagator goes on shell. 
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Figure A. 3: Diagram for m6 — ^ dt{-^ bW' 
For simplicity we have omitted a common factor of 



9s \ 



g^V ^ {2Eu) {2E,) ) {2E,) (2E„) 
V2) {t^-m^){wl-M^){wi-M^) 



(A.39) 



and color matrices.^ We note that the polarization vectors for spin-1 gauge bosons 
may be expressed in terms of spin-| bras and kets. We define 



then for the transverse polarizations 

= TV2\g+){g-\ 
-- ±V2\g-){g+\ 



(A.40) 



(R) 



and for massive spin-1 gauge bosons 



(A.4i: 





m 


E 




\ m 


m 




m 

\ m 


E 
m 



2-\g+){g- 
m 

m 



(A.42) 



A. 4 Helicity Amplitudes for uh ^ dt 

t Vl)) 



In Figure [A.3| , we show the diagram for ub ^ dt{-^ bW'^{-^ i^Ui)). The matrix 
element in the helicity amplitude formalism for this process is 



M = A{u + \bi-){e+\ \ d+){b2-\h+) 



(A.43) 



The color factor for the amphtude squared is3x4x l/3x 1/8 for this process. 
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where we ignore the b parton mass. We have omitted a common factor of 



^V2) {t^-ml){wl-Ml,){wl-Ml,) ^ ■ ^ 



and color matrices.^ 



A. 5 Helicity Amplitudes for ud ^ W* bt{- 



In Figure A. 4 , we show the diagram for ud ^ W* ht{-^ hW^{-^ ^e))- Aside 



from a possible phase, the matrix element in the helicity amplitude formalism for this 
process (W* production) is identical to that of Figure [A.3| . This is because one process 
is the cross diagram of the other and therefore only the momentum assignments are 
different. For clarity, it is 

M =A{u+\bi-){e + \ f^\d+){b2-\n+). (A.45) 
Again we ignore the b parton mass and omit the common factor of 



'g^y ^{2E^){2E,) {2E,, ) {2E,, ){2E,) (2E„) 



(A.46) 



and color matrices.^ 



A. 6 Helicity Amplitudes for ud ^ hhW^{-^ Ui] 



In Figure |A.5| we show the diagrams for ud —>■ bbW^{^ the major W + 2 



jets background to VT-gluon fusion including 6-tagging. The matrix element in the 

^ The color factor for the amphtude squared is 3x3x1/3x1/3 for this process. 
^ The color factor for the amplitude squared is3x3x l/3x 1/3 for this process. 
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Figure A. 5: Diagrams for ud ^ bbW^{- 



helicity amplitude formalism is 



M{+,-) =4:{d-\n+){bi+\u-){e+\ ^-\b2+)/q^ 
+ 4:{d -\b2 +) {e +\ ii -) {bi +\ /- \ n+)/r^ 

Mi-,+) = A{d-\n+){b2+\u-){e+\ 4- \bi+)/q^ 
+ A{d-\bi+){e+\u-){b2+\ \ h+)/r'^ 



(A.47) 



where we have indicated the helicity states of b and b as A1(Ab^, Afej)- We have once 
again left out the factor of 



9s\ 



9w 



{2Eu) {2Ed) {2Eb, ) {2E,,){2Ee) (2E„) 



(A.4J 



and color matrices. [] 



The color factor for the amphtude squared is 2 x 1/3 x 1/3 for this process. 



Appendix B 

Event Rate of the (2 3) Process 

ug dtb 



Monte Carlo integration is an indispensable tool in phenomenology. However, when 
performing a calculation one often encounters singularities which make it impossible 
to obtain meaningful results from a Monte Carlo program. In the case of delta 
functions one is forced to integrate by hand. Other singularities may occur when 
propagators go on mass-shell. These types of divergences may be regularized by 
applying suitable cuts on the external particles in the process. However, when one 
is interested in obtaining a total rate, part of the calculation must be performed by 
hand if there is any singularity present. 

To obtain the total rate for the (2 3) process ug —>■ dtb via Monte Carlo 
integration, the idea is to integrate out the W^g th sub-cross section in Figure ( |B.1| ) 
by hand. The alternative is to use Monte Carlo for the full 3-body phase space, but 
the singularity due to the small mass of the h quark in the right diagram would 
require too much computer time even for the powerful technique applied in VEGAS, 
a fortran code for calculating multiple dimensional integrations [Q. An added benefit 
of performing the sub-cross section integration in the way shown in this appendix is 
an understanding of the vahdity of the effective-VF approximation [R3 . 




Figure B.l: Diagrams for ug — > dth . 
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The differential cross section for this process is 



1 f d^Ti 



where 

ft 

(27r)3(2£;,) (27r)3(2£;,) 



Since_(j5ub is a Lorentz invariant, we choose to evaluate it and the matrix element in 
the tb center of mass (CMS) frame. We express the 4-vector components of t, b, e, n 
and k explicitly as 

t — {Et,p sin 9 cos (l),p sin 9 sin (j),p cos 9) 

b — {Eb, —p sin 9 cos (j), —p sin 9 sin 0, —p cos 9) 

e = sin ^e,0, cos ^e) 

n = sin 6I„,0, cos 6I„) 

k = 0,0,1). (B.3) 

In the tb rest frame s = [t + hY, so 



E. 



s + ml — ml 



2V§ 



s — ml + ml 



p = 



2v^ 

^\l[s - {mt + mb)2][s - {mt - mb)^]. (B.4) 



2v^ 



— — * 

It is easy to see in the tb CMS, that e, k, n all lie in the same plane. Momentum 

conservation {e + k = n + t + b) and t + b = imply e + k = n. It only takes two 
vectors to describe a plane and any vector which is a linear combination of those two 

vectors lies in that plane. Therefore, we choose the vectors e, k, n to define the x-z 

plane with the momentum of the gluon (A;) along the 2;-axis. 

For a 3-body final state, there are only 5 (= 3x3 — 4) independent variables. (The 
minus 4 is for energy- momentum conservation.) We shall choose the 5 independent 
variables to be 9^ (p, En, 9n and 0„ where 0„ can be trivially integrated as / d(f)n = 27r. 
This is a consequence of the arbitrariness of the choice of scattering planes. Hence, 
we should be able to express all vectors in terms of the 4 variables 9, 0, E^ and 9n- 

What we have done so far is express all 4- vectors in the tb CMS. It is now a fairly 
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trivial exersize to evaluate cXsuft.Q We found 



a.„fc = — Id cos ed<f)\M\^. (B.5) 



Having done so, we now have a Lorentz invariant expression for 

^sub ^subi^^j -^n; ^n) ■ 

To compare our result with the effective-W^ approximation, we express En and 6n in 
terms of ek CMS quantities with the aid of some projection operators defined below. 

To extract out the information of e. A; and n in the tb CMS frame we use the fact 

that e + k — n = and {e + k — n)"^ = s. We then define the projection operators P^; 
and P|| as follows: 

P,V . (l±i_!!LZ (B.6) 

Vs 



Ek 

where V = e,n or k. Pe and Py project out the energy of V and component of V 
parallel to k {i.e., the z-axis) respectively. We express the 4- vector components of 
e, n and k explicitly as 

k = ^(1,0,0,1) 
e = ^(1,0,0,-1) 

n = <(l,-sin^;„0,-cos^;,) (B.8) 

where s = {e + k)"^. From now on we will distinguish ek CMS components from tb 
CMS components by a prime. We assume 

< = <(1-^) = ^(1-^) 

and use the following results 

s = (e + k — nY = xs 

= {e-nf = --{l-x){l-cose'j. (B.9) 



^ This integration is done analytically to avoid bad convergence in the numerical integration 
method. All the singular terms are of the nature of In(m^). Because the analytic form for this result 
is long, it will not be explicitly given here. 
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Using the projection operators it is easy to show 



En = —^S{l-X) 

2Jxs 



e|| = E,cosee=^{W^ + s{l-x)) + Er, 

nil = E„cos^„ = -^ + ^e (B.IO) 



and therefore 



Also, the Lorentz invariant phase space integral in terms of ek CMS components 
becomes 

i {27r)%2E'J (47r)2yo ^-.(i-x) ^ ' 

Finally, the following cross section 

1 /-o 



a — 



2s(47r)2Jo i-s(i-x) 



may be folded in with parton distributions and safely evaluated using a Monte Carlo 
program. 



We note that yW^ is the virtuality of the PF-boson line emitted from the w-quark 
line (with momentum e). {1—x) is the fraction of the incoming w-quark energy carried 
away by the outgoing d-quark line in the ek CMS. If desired, one can approximate 
the above equation as the result of the effective-H^ approximation. However, we shall 
not pursue it further here. 



Appendix C 



Helicity Amplitudes of t ^ W^b 

and t W~b 

In Equations ( |6.1|) and (|6.2| ) we have listed the most general form factors for the decay 
processes t + b and i W~ + b. Here we use those equations to calculate 

the helicity amplitudes for an on-shell ly-boson. (We take the limit of — > in 
the following for ^ m^.) 

For the decay process t —>■ W^b, the top quark is taken to decay in its rest frame 
where the top quark momentum is pt = (mt, 0,0,0). Spherical coordinates are used 
to describe the outgoing particles; 6 is taken from the positive Z-axis and (j) is taken 
from the positive X-axis in the X — Y plane. The bottom quark and the W^-boson 
are taken on their mass shells with the four-momenta for the bottom quark {pb) and 
the VT-boson (pw) taken as 



Pb = {Eb, —Eb sin 6 cos (p, —Eb sin 6 sin 0, —Eb cos 6) , 
Pw = {EiY, Eb sin 9 cos (p, Eb sin 9 sin (j),Eb cos 9), (C.l) 

where we have neglected the bottom quark mass, and 

E,^'-^^. (C,2) 

The angles 9 and (p refer to the direction of the W^-boson. 

Denote the helicity amplitudes as {ht, Xw, hb) with A^y = —,+,0 being a left- 
handed, right-handed, and longitudinal ly-boson. After suppressing the common 
factor 



-^^/2Ebmt, (C.3) 

there are 8 nonvanishing helicity amplitudes in the rest frame of the top quark for 
rrib = 0: 
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+ 



- + + 



+ + 



+ + + 



e_ 

2' 



\Mw 



e 
-( 

2 



\Mw 



Ml 



= -V2(/; 



-2i</> 



' fn COS -e 

V 2 



(C.4) 



To obtain the averaged amplitude squared, a spin factor | should be included. We 
note that there is no right-handed VF-boson produced with a massless left-handed 
b from a top quark decay. Similarly, from helicity conservation, it is not possible to 
have a left-handed IV-boson produced with a massless right-handed b from t decay. 

For an unpolarized top quark decay, after summing over the helicities of the bot- 
tom quark, the amplitudes squared for various W polarizations are, apart from a 
common factor (g'^Ehmt), 



\M{Xw - 



'J2 



Ml 



w 



'J'. 



Ml 



w 



Ml 



w 



-ft + 



+ 



Ml 



w 



(C.5) 



The fraction (/Long) of longitudinally polarized VF-boson produced in the rest frame 
of the top quark is defined as the ratio of the number of longitudinally polarized 
VF-bosons produced with respect to the total number of VT-bosons produced in top 
quark decays: 
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r ^ ^i^W = 0) 



|M(Aiy = 0)|2 + \M{Xw = -W+ \MiXw = +W 

where we use T{Xw) to refer to the decay rate for a top quark to decay into a VT-boson 
with polarization Avk- 

Using a parallel definition for the process i —>■ W~b, we obtain the helicity am- 
plitudes {ht, Xw, hi), similar to the ones listed in Equation (|C.4|) for the t — > W^b 

process, provided we replace f[ by //^*, by f[*, /|' by /^*, and by /|'*. (Here 
the superscript * means complex conjugate.) 

The helicity amplitudes of the process e^Ve are well known. After sup- 

pressing the common factor {gMw), they are 



(A^ = -) = -e-^-f ^T'- 



{Xw = 0) 



sing: 



where and 0^ refer to in the rest frame of . 

The helicity amplitudes (A^^) for the decay process W" e"z7g can be obtained 
from Equation (|C.7| ) by replacing by tt — 9* and (f)*^ by tt + </>:. In this case, 9* and 
6* refer to e~ in the rest frame of W~ . 



Appendix D 

The Total Rate for gluon Fusion 



As discussed in Section 2, the total rate for the P^-gluon fusion process is obtained 
by 

Total = (2 ^ 2) + (2 ^ 3) - (sphtting piece) 

and the rates of 

(2^2)=! rf^i d^2 /g'M(6, /i) fb/B{^2, /i) H<l'b ^ qt) 

+ J d^i d^2 fb/Ai^i, IJ) fq'/Bi^2, lA crihq qt) (D.l) 

(2 ^ 3) = y" dil d^2 fg'/Ai^i, /i) fg/Bi^2, /i) ^{q'g ^ qtb) 

+ J d^i d^2 fg/A{.iu /^) fq'/Bi^2, /i) ^{gq qtb) (D.2) 

(sphtting piece) = J d^i d^2 fq'/A{^i, /u) Z^) ^{g'b qt) 

+ J d^i d^2 fb/Ai^i, /i) yu) gt) (D.3) 

where, for instance, /^/^(G,/^) denotes the parton distribution function (PDF) of the 
b quark inside hadron A, carrying the fraction of the hadron momentum, and fi is 
the energy scale at which the PDF is evaluated. The constituent cross section a is 
given by the differential cross section for 

pi(mi) +P2{m2) ^ p'i{m[) + P2{m'2) + ■ ■ ■ +p^(m^), (D.4) 

as shown in Figure |D.1| . The differential cross section is 



da = l-^NM^i+^^^K,---,K) 5) 
4\/(pi ■P2)'^ - mlml 
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with 



i=l 



(D.6) 



and \Ai P is the square of the amphtude after summing over the final state polariza- 
tion and colors and averaging over the initial state polarization and colors. Notice 
that the differential cross section for da{bq' qt) can be obtained from da{q'b qt) 
by interchanging the 4-momenta pi and p2 in the scattering amplitudes. In terms of 
the polar angle 6* and azimuthal angle 0* defined in the center-of-mass frame of the 
partons, this means replace by vr — 6* and 0* by vr + 0*. In Equation ( |D.3D , the 
effective parton density 



(D.7) 



in the MS scheme. The coupling constant 



127r 



(33-2n^)ln(^; 



(D.8) 



Here, nj is the number of quarks with mass less than the energy scale fi. The QCD 
parameter A = A^'' is an experimentally determined parameter. Perturbative QCD 



is presumed to be applicable for A. For CTEQ2L PDF, A- 



MS 



190 MeV. 




Figure D.l: n-body scattering 



Appendix E 



The Eikonal Approximation for 

a{AB{ud) bbW^ + jet) 



Applying the Eikonal approximation, we can write the cross section of a{AB{ud) —>■ 
bbW'^ + jet) in terms of the amplitude square of the ud bbW'^ process as follows. 



a{AB{ud) bbW+ + jet) = JdQ^ JdyJ dql J rf$3 — 



1 f^siOY 

27rq^ 



X 



fu/A{xA,Q) 



{xb,Q) + 



(E.l) 



fu/A{xA,Q) fd/A{xB,Q) 



\M{ud^bbW+) \ + 



{A^B), 

where d^s is the usual 3-dimensional phase space volume as defined in Appendix B. 
Q, y and qt are the invariant mass, rapidity and transverse momentum of the 



_ _ _ 2 — — 

{W + b + b) system. \M.{ud — > bbW+) \ is the amplitude square of ud — > bbW'^ after 
summing over the spin and color factors in the final state and averaging over the spin 
and color factors in the initial state. For a given Q, y and qt 



_ Q y _ Q -y 



where \/S is the center-of-mass energy of the hadrons A and B. 
The splitting functions are 



[E.2) 



(E.3) 
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(E.4) 

and 

(E.6) 

where the indices j and k denote the flavor of quark or antiquark, 6kj is equal to 1 for 

k = j and zero otherwise. In QCD, Cp = | for three colors. The "+" prescription is 
defined by 

where 

^-^ I otherwise ' (^■^) 

The above result holds in the soft-gluon approximation. We have also assumed 
that the initial state QCD radiation dominates the soft-gluon radiation from ud 
bbW~^ . This should be a good approximation because the b quark is massive and is 
less likely to radiate gluons as compared to the initial state quark or gluon. 



Appendix F 

The Computer Program ONETOP 



Our analysis is based on our Monte Carlo program ONETOP, created by modifying 
PAPAGENO (version 3.07), written by Ian Hincliliffe. ONETOP contains code for 
parton level analyses of single top-quark production at hadron colliders as well as 
the major background. The top quark decays on-shell to bW^ with branching ratio 
Br = 1. All final state W^s decay on-shell to ez/ with branching ratio Br = |. 

In addition, we implemented QCD ti production with the top quark decaying on- 



shell according to the Effective Lagrangian of Equation (|6.1|) , which includes the most 
general t-b-W couplings. 



Only the CTEQ2 leading order parton distribution is implemented ||27|| . ONETOP 



accepts matrix elements calculated using the helicity amplitude method described in 
Appendix A. Squaring of the matrix elements and sums over spin and color are 
performed numerically. We include a simple histogramming package which allows the 
plotting of one- and two-dimensional differential cross-sections to aid in analyzing 
event topologies. 

We list below the processes included in ONETOP. The structure of these processes 
are fully discussed in Section 7. 



• q'g qt{^ bw+{^ i+u)) 6, 

• q'q ^W* ^ bt{^ bW+{^ e+u)), 

• q'q^bbW+{^ e+u), 

• qq, 99 ^t{^ bW+{^ e+u))t{^bW~{^ i-jy)). 
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